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SUMMARY: This is the final report for grant NAGW 2055, "Studies of Dynamic 
Processes Related to Active Experiments in Space Plasmas", covering research performed 
at the University of Michigan in the period August 1, 1990 to September 30, 1992. The 
grant is a continuation of grant NAGW 1566, which was awarded to the group while 
located at Stanford University. The grant was awarded to study (a) Theoretical and data 
analysis of data from the CHARGE-2 rocket experiment (IkeV; 1-46 mA electron beam 
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1. Introduction 


The grant NAGW 2055 substituted NAGW 1566 when the PI and his collaborators 
moved from Stanford University to the University of Michigan. The period covered in this 
report correspond to the time at the University of Michigan namely August 1, 1990 to 
September 30, 1992. This period includes a one year no cost extension. 

The work under this grant "Studies of Dynamic Processes Related to Active 
Experiments in Space Plasmas" has covered modeling and data analysis of electron beam 
ejections form the CHARGE-2 sounding rocket experiment (1 keV, 1-46 mA) and the 
Spacelab-2 experiment (1 keV, 100 mA), studies of the dynamics and radiation from 
plasma clouds ejected into ambient magnetized plasmas, and studies of radiation from 
tethered spacecraft systems. At the time of the move to Michigan, the tethered studies were 
completed along with the data analysis of CHARGE-2 data. Also completed was a two- 
stream model for the interaction of electron beams with the Earth’s neutral atmosphere. 

During our time at the University of Michigan, the two-stream model was applied 
to a range of situations, studying the effects of secondary electron fluxes on spacecraft 
charging, density enhancements and optical emissions in the beam column, and electron 
energy spectral variations with altitude. To extend the beam studies, a 3-D electromagnetic 
and relativistic particle code was developed. With the code we have simulated the dynamics 
and the radiation of an electron beam, for the first time in 3-D with a fully electromagnetic 
code (to our knowledge). In addition, we embarked on a study of plasma clouds 
propagating across a magnetic field, inspired by CRIT and CRRES results. More details on 
our work is given in the sections below. 

2. The Beam Atmosphere Interaction 

The Beam Atmosphere Interaction (BAI) has been studied by means of a two- 
stream code. The code solves two first-order, non-linear differential equations in the 
differential energy flux of particles streaming along a magnetic field. The equations couple 
the forward and backward fluxes, and model the transport, secondary electron generation, 
and the cascading in energy of primary and secondary electrons. The code has been 
successful in modeling the return current contribution to the CHARGE-2 rocket resulting 
from secondary and back scattered primary beam electrons. CHARGE-2 is the first 
experiment that has allowed an accurate measurement of the BAI return current 
contributions. The good agreement between model and observations indicate that we have a 
good understanding of the BAI process. 



The model also produces optical emission rates, plasma density enhancements, and 
the energy spectrum of electrons inside the beam column as a function of altitude. This has 
proved very useful in the SEP AC experiment flown on the ATLAS- 1 mission. One of the 
objectives was to generate artificial aurora by a beam ejected from the shuttle and observed 
by cameras in the shuttle bay. The optical emissions observed in the upper atmosphere as a 
result of beam ejections are currently being modeled with the code. 

Our BAI results as well as other advances in the field over the past decade was 
published in a comprehensive review paper [Neubert and Banks, 1992]. It is included in 
the appendix. 

3. The 3-D Particle Code 

To further study the electromagnetic radiation and the dynamics of electron beams, 
a 3-D electromagnetic and relativistic particle code was developed. As a new feature "quiet 
start" conditions were approached by including the poloidal magnetic field due to the 
current carried by the beam electrons streaming along a background magnetic field. The 
complex structure of the whistler mode waves in the vicinity of the beam suggest that the 
transverse velocity (gyration) of the beam electrons, generated as a result of Langmuir 
space charge oscillations, is involved in the generation of the whistler mode. A paper 
describing these initial results has been submitted to Geophysical Research Letters 
[Nishikawa et al, 1992]. The manuscript is included in the appendix. 

The code has also been applied to the problem of plasma clouds ejected across a 
magnetic field. The motivation was here the recent CRRES and CRTT releases. With the 
full 3-D configuration, it was possible to study the formation of charge layers at the sides 
of the cloud, and on the front and rear surfaces of the cloud, and the expansion of these 
layers along the magnetic field into the ambient plasma. It was found that the electrons 
accelerated out of the charge layers may reach very high energies and directly generate 
optical emission in the upper atmosphere (aurora). Results on the dynamics of plasma 
clouds and the electromagnetic radiation from such clouds was published in Journal of 
Geophysical Research [ Neubert et al., 1992]. The paper is included in the appendix. 

Finally the code was used to simulate a plasma flowing over a magnetic dipole, 
analogous to the solar wind flowing past the Earth's magnetic field. The results show the 
formation of a bow shock an extended tail, the Chapman-Ferraro current systems, and the 
penetration of solar wind particles into the cusp regions. The paper will be published in 
IEEE transactions on Plasma Science [Bunerman et al., 1993]. It is included in the 
appendix. 
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RECENT RESULTS FROM STUDIES OF ELECTRON 
BEAM PHENOMENA IN SPACE PLASMAS 

TORSTEN NEL BERT and PETER M. BANKS 

Space Physics Research Laboratory, Department of Oceanic, Atmospheric and Space Science, 
University of Michigan, Ann Arbor, MI 48109-2143, USA. 

( Received 1 7 July 1991) 

Abstract — Experiments involving the ejection of beams of electrons from spacecraft have been performed 
for more than 2 decades in order to study fundamental plasma physical processes as well as for a range of 
diagnostic- and application-oriented purposes. This paper reviews some of the key issues that have been 
pursued in the past 10 years. These include questions regarding spacecraft charging and beam dynamics, 
the interaction of beams with neutral gas and plasmas, and the electromagnetic radiation generated by 
continuous and pulsed electron beams. It is shown how our understanding of these phenomena has 
matured, thereby providing a solid foundation for future experiments involving the use of electron beams. 


1. INTRODUCTION 

The Universe as we know it is made up of neutral 
and ionized matter permeated by gravitational and 
electromagnetic fields. It is highly dynamic and in con- 
stant change as manifested by the birth and death of 
stars. One of the fundamental processes that arises 
naturally out of this dynamic interaction between 
matter and fields is the generation of electric currents. 
As charged particles tend to flow in force-free direc- 
tions, currents tend to be magnetic field-aligned in 
which case they are often called “Birkeland" currents 
after Professor Birkeland who first suggested such 
currents flow in the Earth's magnetosphere (Egeland 
and Leer, 1986). Their existence is now firmly estab- 
lished and it has been found that they are intimately 
connected to “inverted V" events, the generation of 
aurora, and the generation of electromagnetic radi- 
ation in the range from a few r Hertz to many mega- 
hertz. It is also thought that Birkeland currents are of 
importance throughout the Universe on a variety of 
scales up to galactic scales (Faelthammar, 1986; 
Peratt, 1986, 1988, 1990). It is therefore not surprising 
that the scientific community has shown considerable 
interest in performing space experiments which 
involve the ejection of electron beams and the for- 
mation of artificially generated currents. These so- 
called “active experiments" allow the study of fun- 
damental plasma physical processes, thereby gaining 
quantitative information that is almost impossible to 
obtain when observing natural phenomena. With 
active experiments, the charged particle source (e.g. 
the electron beam) characteristics are known and pro- 
vide a specific stimulus to the coupled response of 
electromagnetic fields and plasma. In many ways. 


active experiments provide an important calibration 
of natural phenomena. 

There are a variety of additional objectives that 
also motivated past space experiments with charged 
particle beams. These include the probing of the 
Earth's magnetic field on a global scale by recording 
the echo of particles mirroring in the conjugate hemi- 
sphere as done in the ECHO series of sounding rockets 
(Winckler et ai , 1989; Hallman et ai , 1990), probing 
the local electric and magnetic field as done on GEOS 
1 and 2 (Melzner *>/<;/., 1978 ; Junginger et ai, 1984), 
and using pulsed electron beams as a means of gen- 
erating electromagnetic radiation for purposes of 
communication between ionospheric spacecraft and 
receivers on conjugate spacecraft or on the Earth (e.g. 
Inan et ai , 1984). 

However, experience has shown that there are many 
technical problems associated with the successful ejec- 
tion of charged particle beams in space and the study 
of these has taken considerable effort. A foremost 
challenge has been the charging of spacecraft plat- 
forms as a result of charged particle beam ejection. 
As pointed out by Linson (1982). an isolated spherical 
conductor of 1 m diameter will charge at a rate of 9 V 
/is 1 when a beam of 1 mA is ejected and no return 
currents from the ambient plasma are collected. 
Because of this rapid charging rate, and its influence 
on processes in the plasma surrounding the beam 
source, it is important to undertake experiments and 
theoretical studies aimed at understanding and miti- 
gating the processes associated with electrical charg- 
ing of the beam-emitting platform. In doing so the 
goal is both scientific and technical. On the one hand, 
to design successful space experiments with a charged 
particle beam apparatus it is important to understand 
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the underlying physical processes linked with the 
beam emissions. On the other hand, it is also impor- 
tant to be able to achieve a high quality of the emitted 
charged particle beam for the purpose of probing 
nature with a beam of well-defined characteristics. 
The quest for such understanding involves an exciting 
body of basic physical processes, including the for- 
mation of large and complicated sheath regions of 
charged particle and electromagnetic field fluc- 
tuations, plasma-neutral gas interactions, wave- 
particle interactions within plasmas, and non-linear 
phenomena such as plasma wave turbulence, just to 
name a few . 

In this paper it is our aim to review selected results 
of the 1980s and to discuss future directions of active 
electron beam experiments in space. More infor- 
mation on beam experiments is found in the review 
paper by Winckler ( 1 980), the monograph by Grandal 
(1982), the special issue of Advances in Space 
Research , edited by Brenning and Mendillo (1990), 
the special issue of Radio Science (1984). and the 
European Space Agency report by Burke (1983). 

The organization of this paper is as follows : in 
Section 2 some basic processes associated with 
electron beam emissions are discussed. These are 
spacecraft charging, including the phenomenon of 
super charging and virtual cathode formation, beam- 
atmosphere interaction (BAI). and beam system- 
plasma interaction (BSPI). Then, in Section 3 we dis- 
cuss the theory for electromagnetic radiation from 
continuous and pulsed electron beams and relate some 
experimental results obtained from space shuttle and 
sounding rocket experiments. Finally in Section 4 we 
discuss the current status of electron beam exper- 
iments and point out future directions for the field. 

2. BASIC PROCESSES 

When an electron beam is ejected from a spacecraft, 
the spacecraft potential will rise, drawing return cur- 
rent electrons from the ambient plasma. For low 
potentials, space charge in the plasma is small and the 
electrons collected by the spacecraft are those whose 
trajectories directly intersect the spacecraft structure. 
For potentials well above the ambient plasma thermal 
energy, an extended region of positive space charge 
will develop around the spacecraft. Ambient electrons 
that enter this sheath have a chance of reaching the 
payload through the combined action of the highly 
inhomogeneous electric field in the sheath and the 
ambient magnetic field. Thus, the effective current 
collecting area of the spacecraft is enhanced by the 
space charge sheath, thereby increasing the collected 
current. 


The electron thermal current density / iK in an equi- 
librium plasma is given by : 

Jo< = e«<r c > 4. (I) 

where e is the elementary charge, n the plasma density, 
and <r t .) the average electron velocity. In a Maxwell- 
ian plasma the expression for <r c > is : 

<r c > = (8kr/Jim c )' 2 , (2) 

where kT is the electron energy and m e the electron 
mass. The ion thermal current density and ion velocity 
are given by similar expressions. While the tem- 
perature of the ionospheric plasma is fairly constant 
and approximately the same for ions and electrons to 
within a factor of 2-3, the densities may vary many 
orders of magnitude. Typical values for the particle 
velocities and the peak densities encountered in the 
ionosphere and the corresponding current densities 
are given in Table 1 . 

The electron thermal current available to a sound- 
ing rocket payload at low potentials in the ionosphere 
is 2-10 mA 2 depending on the plasma density and 
temperature. Beam experiments flown in the past have 
carried beam accelerators ejecting beams with cur- 
rents up to several hundred milhamperes and indi- 
vidual electron energies up to 40 keV. Since current 
collecting areas on such platforms are seldom more 
than a few square meters, such payloads have often 
reached high positive potentials. 

For most applications it is essential that the beam 
escapes the nearby region of the spacecraft in a well- 
known and coherent manner, preferably with the full 
electron accelerator energy and small energy, spatial 
or temporal dispersion. In practice, observations and 
theoretical simulations indicate that beam emissions 
are accompanied by many dynamic charging 
processes, both of the platform, w ithin the beam itself, 
and in the surrounding medium. Thus, platform 
charging and beam characteristics can have large 
dynamic fluctuations on time-scales of the order of 
plasma time periods (microseconds) and physical 
scales of the order of Debye lengths (millimeters). 
These micro-scale effects have the net consequence of 
d.c. charging of the spacecraft, reducing the beam 
escape current, degrading the average beam particle 
energy by an amount corresponding to the average 
spacecraft potential, and degrading the beam coher- 
ence. 

A number of effects help to enhance the emitting 
platform return current above the limit set by iono- 
spheric thermal current densities and thereby reduce 
platform charging. These effects are associated with 
interactions between the ambient plasma and any neu- 
tral gases surrounding the payload and the electron 
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beam and charged platform. Specific effects include 
beam-atmosphere interactions (BAI) and beam-sys- 
tem plasma interactions (BSPI). The neutral gases 
involved may be the ambient atmosphere, outgassing 
from the platform, or attitude control thruster emis- 
sions. By various processes, these gases can be ionized 
and thereby create enhanced plasma densities in the 
vicinitv of the payload. This component of the return 
current is sometimes called the “active component 
to distinguish it from the “passive" component which 
is the current drawn to an anode immersed in a 
plasma. The process identified as the beam-plasma 
discharge (BPD). investigated in many chamber 
experiments (Bernstein el a/., 1979), is thought to 
occur in space at relatively low altitudes (110-160 
km) (Linson. 1982) where the neutral gas density is 
relatively high. It will not be covered in this review. 

2.1 . Passive current collection 

The problem of determining the current and poten- 
tial existing between two concentric electrodes in an 
unmagnetized plasma was first discussed by Langmuir 
and Blodgett (1924). In their theory, an anode 
immersed in a plasma will develop a positive space- 
charge sheath out to a distance r s , the sheath edge, 
where the potential reaches the ambient plasma poten- 
tial. At the sheath edge, thermal electron fluxes enter 
the sheath from the ambient undisturbed plasma, and 
the potential structure in the sheath as well as the 
location of the sheath edge are calculated by con- 
sidering current continuity and Poisson's equation. 
The mixed boundary conditions are the thermal cur- 
rent density at the sheath edge and the potential of 
the conducting bodies. From this theory the current 
to a probe is space -charge limited. Langmuir and 
Blodgett (1924) gave the expression for the current to 
a probe at a positive potential <J> 0 in the case of spheri- 
cal symmetry and no magnetic field : 

4, = 47i£ 0 (32e/8 lm c ) 1 2 <D^ "a" (3) 

where a is a geometrical parameter. It was assumed 
that the current was carried by electrons to an anode. 
If the probe sheath radius is much larger than the 
probe radius. r p , the following expression was found 
for 2 : 

a 3;4 « l.llr,/r p . (4) 

Since the current to the probe must equal the cur- 


rent across the probe sheath, we also have : 

/*! = 4tt rljoc. (5) 

Combining equations (3)— (5) we find : 

/„, = CjL\r p < D 0 ) 6 7 . (6) 

The expression for the ion current to a cathode is 
similar. For electrons the numerical value for C c is 
6.946 x 10 \ and for 0 + ions C, is 3.674 x 10 ~ 4 , with 
the parameters of the equation given in SI units. 

The current given in equation (6) has been found 
for an unmagnetized plasma. This condition can be 
regarded as a rather good approximation when con- 
sidering ion collection by a negatively biased payload 
because of the large ion Larmor radius. For example, 
an 0 + ion which has just entered the sheath and 
obtained the energy of 1 eV will already at that point 
have a Larmor radius of 13m (B = 0.44 G) which 
usually is large compared with the sheath dimension. 
The ions can therefore be regarded as unmagnetized 
in the sheath region. Electrons are, however, limited 
by the magnetic field in their motion perpendicular to 
the field, and thus the electron current that can be 
collected is reduced compared with the values pre- 
dicted by equation (6). 

Theories extending the work of Langmuir and 
Blodgett (1924) for the case of beam-emitting space- 
craft in the ionosphere was put forward by Beard 
and Johnson (1961) (without magnetic field) and 
Parker and Murphy (1967) (with magnetic field). 
Parker and Murphy included the magnetic field effects 
by using the conservation of energy and angular 
momentum of particles collected by a sphere. The 
upper limit of electron current collection according to 
their results is: 

/// 0 < 1 + [8e<D 0 i(m c a)ia 2 )] 1 2 , (7) 

where 7 0 is the thermal electron current to the sphere 
at plasma potential, is the electron gyrofrequency 
and a is the radius of the sphere. 

The space-charge limited current for an unmag- 
netized plasma given in equation (6) increases as O 6 
while the magnetic-field limited current in equation 
(7) increases <J>‘ 2 . A consequence of this and of the 
differing dependence of system parameters appearing 
in the two expressions, the predicted current collection 
in the presence of a magnetic field is smaller than 
without a field as expected. 
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Studies have been made in the past decade in space 
and in laboratory plasma chambers to test the validity 
of the estimates given in equations (6) and (7). For 
space experiments this task is not as simple as one 
would imagine because of the difficulty in measuring 
the spacecraft potential accurately during beam emis- 
sions. As we will show later, the ejection of charged 
particle beams has the effect of generating a turbulent 
region with enhanced energetic particle fluxes in the 
vicinity of the spacecraft. These fluxes contaminate 
the potential measurements, which are usually per- 
formed by measuring the potential difference between 
the rocket body and smaller spheres mounted on 
booms extending from the rocket. Furthermore, in 
many cases the sheaths are of such dimensions that 
the length of the floating probe booms is smaller than 
the sheath radius and thus the potential measurement 
does not reflect the complete potential difference 
between the ambient plasma and the spacecraft. 
Recently however, a technique has been developed 
which uses two payloads connected by a conducting, 
insulated wire. The payloads are allowed to drift apart 
up to many hundreds of meters and the potential of 
the beam-emitting payload is measured with respect 
to the remote tethered payload located outside of the 
disturbed region. 

Such a technique was used in the CHARGE series 
of sounding rocket experiments (Sasaki el al ., 1987, 
1988; Myers et al 1990). The sheath formation and 
the ion current collection to a negatively charged pay- 
load was studied by applying a potential difference of 
up to 460 V between the two tethered payloads, the 
mother payload biased negative and the daughter pay- 
load biased positive. The potential of each payload 
relative to the ambient plasma was estimated by not- 
ing that the total current to the structure must equal 
zero. As a consequence, the ion current to the negative 
mother balanced the electron current to the positive 
daughter. Since the thermal ion current density in 
the plasma is much smaller than the electron current 
density (Table 1 ) the mother charged to a large nega- 
tive potential while the daughter charged to a small 
positive potential. In other words, by far the most of 
the 460 V potential drop was found in the positive ion 
sheath surrounding the mother payload. 

The data obtained in the CHARGE - 2 experiment 
was compared with computer calculations by the 
NASCAP/LEO code (Mandell et al ., 1990; Neubert 
et al 1990a). The code, in essence, simulated the 
space-charge limited current and the corresponding 
sheath structure as given by the theory of Langmuir 
and Blodgett (1924). However, the advantage of the 
code is that the actual shape of the spacecraft can be 
modelled rather than assuming an equivalent sphere 


or an infinite cylinder, as must be done w hen treating 
the problem analytically. Furthermore, the velocity of 
the spacecraft through the plasma could be included 
in the calculation. This is particularly important for 
ion-current collection since the spacecraft velocity is 
of the order of the ion thermal velocity and therefore 
significant wake and ram asymmetries are created in 
the ambient plasma density and the current collection. 
With the computer simulation it was also possible to 
estimate the effects of secondary electron emission 
caused by the bombardment of the rocket skin by 
incoming ions. 

A contour plot of the sheath structure around the 
mother payload computed by the NASCAP/LEO 
code is shown in Fig. la and b for two values of the 
spacecraft potential: —30 and —200V. The contour 
marked with the number 2 is the — 2 V contour, which 
is basically the sheath edge. The length of the rocket 
is 3.26 m and the diameter is 44 cm. Four floating 
probes separated by 25 cm were mounted on a boom. 
The influence of these was also modelled by taking 
into account the internal resistance of the floating 
probe array. The probes can be seen at z = 9 extending 
outward from the rocket. 

The return current characteristic of the mother pay- 
load as a function of rocket potential was determined 
from the measurement of the current in the tether with 
gradually increased bias potentials applied between 
the mother and the daughter. The results are shown 
in Fig. 2 for two such experiments performed at 244 
and 247 km altitude. Also shown are the results of the 
NASCAP/LEO code which show very good agree- 
ment with observations. 

The study of current collection and sheath for- 
mation of payloads at negative potentials led to the 
following conclusions : 

(1) Magnetic field effects on the ion motion in the 
sheath could be ignored, as expected. The sheath 
structure and the current collection was close to the 
predictions of the theory of space-charge limited flow 
with modifications listed below. 

(2) The velocity of the rocket had the effect of 
increasing the ion-current collection by as much as 
40% at low potentials, decreasing to 15% at higher 
potentials (Mandell et al ., 1990). At the time of the 
experiments the rocket velocity was 580 m s" 

(3) Laboratory experiments on secondary electron 
yields arising from alkali metal ions impacting on 
aluminum oxidized surfaces show that significant 
numbers of secondary electrons are emitted for poten- 
tials above a threshold of 200 V and that the yield 
reaches 0.5 electrons per impacting ion at 400 V (Dietz 
and Sheffield, 1975). To the best of our knowledge. 
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Fig. 1. NASCAP/LEO simulation of the potential structure around the CHARGE-2 mother 
PAYLOAD WHEN BIASED TO (a) - 30 V AND (b) - 200 V. 
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Fig. 2. The electron return current -payload potential 

CHARACTERISTIC OF THE CHARGE- 2 MOTHER PAYLOAD (FROM 
N EGBERT el ai, 1990a). 



Potential (Volts) 

Fig. 3. Current collection by the CHARGE- 2 daughter 

PAYLOAD VS VEHICLE POTENTIAL COMPARED WITH : THE LaNG- 

muir and Blodgett model (solid line) and the Parker 
and Murphy model (dashed line) (from Myers et al. , 
1989). 


no studies have been performed for 0 + ion beams 
bombarding oxidized metal surfaces. However, 
reports of other combinations of ion beams and sur- 
faces indicate that the values given above are appli- 
cable to oxidized metal surfaces (Szapiro et ai , 1988; 
Szapiro and Rocca, 1989; Utterback and Miller, 
1961 ; Hayden and Utterback, 1964 ; Amme, 1969). A 
review is given in Langley et al. (1984). The NASCAP 
LEO simulations and comparison with experimen- 
tal data indicated that secondary electrons were emitted 
during the negative bias operations of CHARGE - 2 
roughly as predicted by laboratory data. 

An electron beam source mounted on the 
CHARGE-2 mother payload emitted 1 keV electrons 
with currents up to 46 mA as measured by on-board 
diagnostics. The current-voltage characteristic of the 
daughter payload (electron collection at positive 
potential) was obtained for beam-emission sequences 
executed with low resistance in the tether circuit. Com- 
parison of the observations and the theoretical esti- 
mates given by equations (6) and (7) was reported by 
Myers et al. (1989). The key results are shown in Fig. 
3. The model estimates were found by approximating 
the short cylindrical daughter payload (44 cm diam- 
eter, 82 cm long) with an equivalent (equal area) 
sphere. As can be seen from the figure, the observed 
electron return current and the Parker and Murphy 
model estimates arc in close agreement and sub- 
stantially smaller than the Langmuir-Blodgctt pre- 
dictions. 

With the assumption of current continuity (ejected 
current equals collected current) and beam escape, the 


return current-voltage characteristic of the mother 
payload could also be determined. The observations 
were compared with the models of equations (6) and 
(7), again approximating the mother payload as an 
equivalent sphere. It was found that at high altitudes 
(above 240 km) the electron return current was again 
well below the estimate of Langmuir and Blodgett, 
but above the upper limit of Parker and Murphy 
(Myers et ai , 1989, 1990; Mandell et ai , 1990). At 
low altitudes, however, the return current increased 
above even the Langmuir and Blodgett model esti- 
mate. The reason for this is the increased influence of 
the ambient neutral atmosphere and the BAI process, 
to which we will return in the following subsection. 

Mandell et ai (1990) studied the return current 
collection of the CHARGE- 2 mother payload at high 
altitudes using the true shape of this body. These 
computer estimates were in better agreement with the 
observations than estimates based on the Parker and 
Murphy model. Mandell et ai concluded that the 
difference arose from the fact that the Parker and 
Murphy model assumes a spherical body. Thus, 
although the daughter payload was a short cylinder, 
it was well approximated by a sphere. The mother 
payload, in contrast, was a long cylinder which was 
oriented perpendicular to the ambient magnetic field 
and this significantly increased the return current over 
the Parker and Murphy model values. 

Study of current collection during application of 
passive positive potentials led to the following con- 
clusions : 

(4) For spherical payloads, the electron return cur- 
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rent i> magnetically limited and close to the Parker 
and Murphy limit even during simultaneous electron 
beam emissions, provided the experimental altitude is 
above 240 km. (Rocket apogee was at 260 km so no 
information pertaining to the topside ionosphere was 
obtained.) 

(5) A non-sphcrical payload can collect electron 
return currents somew hat in excess of the Parker and 
Murphy limit, as estimated by approximating the pay- 
load with an equivalent area sphere. The difference 
can be attributed to the fact that the Parker and Mur- 
phy model relies on the assumption of symmetry in 
the ambient magnetic held. 

Laboratory and rocket experiments on charging 
phenomena performed in Japan have been reviewed 
by Kawashima (1982). The observations presented 
there are generally in accord with the conclusions 
reached above, provided the system is in a regime 
where surrounding neutral gases are not ionized. As 
will be discussed in the following sections, neutral 
gases may, in fact, play a very important role in 
space-borne electron beam experiments. 

2.2. Beam-atmosphere interactions 

The beam-atmosphere interaction (BA1) involves 
the collisional ionization of the ambient atmosphere 
by energetic electrons. It is a process which occurs 
naturally on auroral held lines where electron beams 
precipitate into the upper atmosphere generating aur- 
ora and enhancing the ambient plasma densities. The 
excitation of aurora is one of the natural phenomena 
that can be studied with the use of artificial electron 
beams, as will be done in the ATLAS-] space shuttle 
mission scheduled for early 1992 (Taylor et ah. 1991). 

For the purpose of studying the effect of BA1 on 
charging levels and return current collection of space 
craft, the CHARGE-2 sounding rocket experiment 
was the hrst to allow separate measurement ol the 
active and the passive component of the return current 
(Neubert et aL 1990b). The experimental con- 
figuration is show n in Fig. 4, which illustrates the two 
payloads and the current system during electron beam 
ejections from the mother. From the measurements 
of the beam current, and the tether current. Aether* 
the return currents to the payloads can be determined 
by assuming that the beam current escapes the near- 
environment of the mother : 

A> = Aether 

” Avam Aether* 

where /,, is the return current to the daughter and 
/ N1 the current to the mother. The daughter current 
represents the passive current component since the 


( 8 ) 

( 9 ) 



Fig. 4. Configuration of the CHARGE-2 payload and 

THE AVERAGE ELECTRON CURRENT SYSTEM AROUND THE PAY- 
LOAD DURING AN ELECTRON BEAM EXPERIMENT (FROM NEUBERT 

eta!.. 1990b). 


daughter is located far from the mother and outside 
the disturbed region surrounding the mother and the 
electron beam, while the current to the mother rep- 
resents the sum of the passive and active components. 
Since the passive component can be deduced from the 
daughter measurement, the active component can be 
extracted from the mother return current. 

The observations of current collection during the 
CHARGE- 2 flight are summarized in Fig. 5. This 
shows the fraction of the beam current collected by 
the daughter as function of altitude. The labels SQ2- 
SQ6 mark the beam-emission sequences performed 


CHARGE -2 



Fig. 5. CHARGE- 2 observations of tfie tether to beam 

CURRENT RATIO AS A FUNCTION OF ALTITUDE. 

Also indicated is the ratio of the daughter area to the total 
area of the two payloads (from Neubert et at., 1990b). 
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during the Bight numbered in lime sequence. Also 
indicated is the h\ed ratio of daughter current col- 
lection area to the total collecting area of the two 
payloads. At high altitudes, the fraction of the return 
current collected by the daughter is close to the ratio 
of the two areas. This is expected for passive current 
collection. At lower altitudes, a progressively smaller 
fraction of the total return current is collected by 
the daughter. Assuming that the beam escapes, this 
implies that progressively more current must be col- 
lected directly by the mother. This can only happen if 
the active component of electron current increases. 

Since the active current component increases with 
decreasing altitude it was suggested by Neubert et 
al. (1990b) that the active component was due to 
ionization created in a BAI process. In order to quan- 
tify the increased current due to this process, a com- 
puter code, first developed for the study of the gen- 
eration and transport of photoelectrons in the upper 
atmosphere (Banks and Nagy, 1970) and later for 
auroral electron fluxes (Banks et al. y 1974), was modi- 
fied by Neubert et al. for the purpose of modelling 
the behavior of magnetic field-dominated electrons 
ejected from a spacecraft and scattered back and forth 
along magnetic field lines by collisions with atmo- 
spheric gases. Primary and secondary electrons scat- 
tered back to the payload would constitute the active 
source of return electron current. 

The code solves two first-order, non-linear 
differential equations in the differential energy flux 
streaming along an ambient magnetic field. The equa- 
tions coupling the forward flux, with the back- 
ward flux. O , are described by : 

<J> * (r. £) = -<y 2 (r. (r,£) + <r i (r,£)4> (r,£) 

-f<r<r,£)-L(£)«^(r,£) (10) 

<t> (r. £) = <x : (r. £)<D <r.£) -<r,(". £)<*>" U, £> 

-Q <r, £) + £(£)<& (r,£). (11) 

Here O' denotes the derivative of O with respect to z, 
w hich is the dimension in the direction of the magnetic 
field. o : is the cross-section describing the loss in the 
electron energy range £ to £ + d£, a, is the cross- 
section describing the elastic back scattering of elec- 
trons. Q is the electron production rate due to ion- 
ization and cascading of electrons down in energy, 
and L describes the loss of flux out of the beam flux- 
tube induced by the motion of the spacecraft per- 
pendicular to the magnetic field. 

The problem is treated in one dimension and is 
applied to the situation of electrons emitted from a 
spacecraft by modelling the fluxes that are confined 
to the magnetic flux-tube encompassing the beam. It 


is assumed that the beam and the secondary electrons 
fill a flux-tube with a certain cross-sectional dimension 
and that the fluxes are homogeneous within this flux- 
tube. From experiments, an estimate of the width of 
the beam flux-tube is 4 beam electron Larmor radii at 
90 pitch angle (Frank et al 1989; Winckler et al ., 
1989). The loss term L describes the loss of flux out 
of the sides of the flux-tube caused by the motion of 
the spacecraft perpendicular to the magnetic field. 

To model typical conditions, equations (10) and 
(11) were solved numerically for a 1 keV electron 
beam of 100 mA current emitted downwards from an 
altitude of 260 km over White Sands, NM, at 
midnight. The MSIS/86 model (Hedin. 1987) was used 
for the neutral atmosphere (N>. 0 : , O) and the IRI 
model (Bilitza, 1986) was used for the ionosphere. 
The IRI model is needed to model Coulomb collisions. 

The differential electron energy flux at the space- 
craft altitude is shown in Fig. 6. The upward flux is 
the flux returning to the spacecraft from below. A 
portion of this flux will in general be collected by 
the spacecraft and form the active component of the 
return current, while the rest will reach locations 
above the spacecraft and interact with the neutral 
constituents here. In Figs 6 and 7 it is assumed that 
the portion of the flux collected by the spacecraft can 
be ignored relative to the total flux. 

The downward flux consist of two contributions. 


Charge 2: Altitude = 260.0 km 



Fig. 6. Model estimates of the upward and downward 

DIFFERENTIAL ELECTRON FLUX AS A FUNCTION OF ENERGY AT 
THE SPACECRAFT ALTITUDE. 

The beam is emitted downward parallel to the magnetic field 
with an energy of 1 keV and a current of 100 mA (the fluxes 
scale linearly with beam current). The primary beam is seen 
as the peak at 1 keV in the downward flux (from Neubert et 
al., 1990b). 
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One is the electron beam flux. which is modelled as 
Gaussian in energy around the l keV beam encrg>. and 
the other is the secondary electrons generated at lower 
energies and propagating in the same direction as the 
beam electrons from locations above the spacecraft. 
Thus, one will in general expect a hot component of 
electron fluxes to impinge on the spacecratt lrom both 
hemispheres. This is in accordance with observations 
made in the ECHO - 2 experiment (Winckler et al.. 
1975). 

Figure 7 shows the fluxes as a function of altitude 
for 1 keV and 100 eV electrons. As can be explained by 
electron free-path considerations, energetic electron 



Fig. 7. Model estimates of the differential electron 

FLUXES AS A FUNCTION OF ALTITUDE FOR THT SAME PARAMETERS 

used in Fig. 6. 

(a) I-keV electrons and (b) 10-eV electrons (from Neubert 
el a L 1990b). 


fluxes are generated along the magnetic field to con* 
siderable distances from the spacecraft. This is in 
accordance with optical observations in other exper- 
iments w hich found streaks of optical emissions aligned 
with the magnetic field and extending far from the 
payload and its sheath (O'Neil et al.. 1978a, b; 
Winckler et al., 1989). 

At low altitudes the flux becomes omnidirectional 
and approaches zero because of the progressively 
shorter mean free path in the lower thermosphere. 
Above the spacecraft, the upward flux decreases rela- 
tively slowly with altitude because of the increased 
mean free path which allows the electrons to escape 
without significant collisions. Similarly, the down- 
ward flux from above the spacecraft increases from 
zero at high altitudes (a boundary condition) to large 
values at low altitudes. 

The active return current estimated to reach the 
CHARGE - 2 rocket mother payload from such model 
calculations has been found to be consistent with the 
observed magnitude of the active current. Thus, there 
was no need to involve other processes for this exper- 
iment. However, at higher altitudes and for spacecraft 
velocities across the magnetic field typical of orbiting 
platforms, the active return current component due 
to BAI becomes small. In order to estimate the BAI 
currents to other types of ionospheric spacecraft a 
series of calculations were performed for a range of 
altitudes, spacecraft sizes, and velocities (Neubert and 
Banks, 1990). Noting that the BAI return current is 
directly proportional to the emitted beam current, the 
BAI effect can be quantified by a gain factor A which 
is the ratio of the return current to the ejected beam 
current at the spacecraft altitude. The loss to the flux- 
tube incurred by the motion of the spacecraft con- 
tained in the term L of equations (10) and (11) can 
be described by the parameter v : 

v = tv dr. (12) 

where v L is the cross-field velocity of the spacecraft 
and d y is the cross-field dimension of the beam dis- 
turbance in the direction of motion. 

A summary of the BAI study is shown in Fig. 8 
where the gain factor A is plotted as function of v for 
three different altitudes. The gain factor is split into 
two parts describing the fluxes from above and below r 
the spacecraft. As expected, when v becomes large, A 
decreases. Similarly, when the altitude increases, A 
decreases. At the bottom of the figure is indicated the 
range in v that can potentially be covered by various 
ionospheric spacecraft. An array of 10 electron guns 
along the space station may, when the long axis of the 
station is aligned with the velocity vector, result in a 
large value of dy, say 100 m. The ionospheric orbital 
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Fig. 8. Model estimates of the beam atmosphere inter- 

A ( TION ( B A I ) GAIN FACTOR A AS A FUNCTION OF THE 
SPACECRAFT PARAMETER V (FROM NELBERT AND BANKS. 1990). 


velocity is typically 7 km s' 1 with perpendicular vel- 
ocities down to 3.5 km s' 1 . The minimum v in this 
case is 35 s l . 

Since the spacecraft potential sheath is generally 
smaller than the beam flux-lube dimension, A must 
be somewhat larger than 1 in order for the BAI return 
current to balance the beam current. From Fig. 8 we 
estimate that at 200 km altitude, a spacecraft with v 
of 50 s 1 will just be able to collect the complete return 
current from the BAI. 

It is also the case that an energetic electron wake 
exists behind an orbiting or otherwise moving payload 
emitting a steady electron beam. This wake is aligned 
along the magnetic flux-tubes traversed by the beam. 
The width of the region is of the order of four gyro- 
radii. based on measurements by Frank et al. (1989) 
and Winckler et al. (1989). The length of the wake 
behind the payload depends on the time taken for the 
atmospherically scattered electrons to lose their excess 
kinetic energy or to escape from the atmosphere. The 
energy degradation time. 7, depends on the mag- 
nitude of the inelastic cross-section, er,. the neutral 
gas density, h, and the electron velocity. r c . according 


to the relation 7^5 W7 ( r c , assuming that about tne 
collisions arc needed to bring an average energetic 
electron down to thermal electron energies. For con- 
ditions al 220 km, for example, 7 0.03 s. For an 

orbiting payload at this altitude in the ionosphere, the 
energetic electron wake zone would extend approxi- 
mately 1.2 km behind the vehicle. At higher altitudes, 
the wake length is determined by the speed at which 
electrons leave the atmosphere. In this case elastic 
scattering dominates, but numerical models are 
necessary to determine counteracting effects of the 
atmospheric albedo for downward moving electrons 
and net escape from the atmosphere. 

The plasma density enhancement and optical emis- 
sions from the BAI are phenomena which will be 
explored in the A TLAS - 1 space shuttle mission (1992) 
and the CHARGE- 2b sounding rocket experiment 
planned for the spring of 1992. The experiments can 
study auroral processes in a quantitative manner and 
can in addition give new r information on the physics 
of plasma density structures in the lower ionosphere 
and their influence on radio wave propagation (Banks 
and Gilchrist, 1985). Using the model described by 
equations (10) and (11), both plasma density enhance- 
ments along the beam flux-tube and optical emission 
rates can be calculated. Examples of the plasma den- 
sity enhancement created by electron beam ejection 
from 200 km altitude are shown in Fig. 9 for two 
electron beam energies. The time constant for recom- 
bination is inversely proportional to the electron den- 
sity and is about 7.7 s for 10 6 cm 3 increasing to about 
1 min when the density has decayed to 10' cm At 
this point the density is still a factor 100 larger than 
the background density. 

In situ observation of plasma density enhancements 
is difficult. However, observing the phase shift of the 
telemetry signal between a mother and daughter pair, 
where one of the payloads carries an electron accel- 
erator, estimates of the local plasma density enhance- 
ment can be made as was done in the MAIM IK exper- 
iment (Friedrich et al 1991). The upper estimates 
found here are consistent with the estimates show n in 
Fig. 9. 

The results from the BAI studies were 

(6) The interaction of artificial electron beams w ith 
the neutral atmosphere can be modelled, thereby giv- 
ing quantitative estimates of enhanced return 
currents, as well as plasma density enhancements, 
optical emissions, electron differential energy fluxes, 
etc. along the magnetic flux-tubes containing the 
source electrons. 

2.3. Beam! system-plasma interactions 

The ejection of an electron beam from a spacecraft 
strongly perturbs the local plasma environment. In 
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Fig. 9 Model estimates of the plasma density in the 

BEAM FLUX -TUBE AS A FUNCTION OF ALTITUDE (HEAVY LINE) 
AND AMBIENT PLASMA DENSITY (THIN LINE). 

The spacecraft altitude is indicated by the dashed line (from 
Neubert and Banks. 1990). 


addition to the ionization of the ambient neutral 
atmosphere mentioned above, beam ejection results 
in plasma wave generation and plasma heating in a 
region surrounding the spacecraft and the beam flux- 
tube (Cartwright et al 1978 ; Amoldy and Winckler, 
1981 ; Arnoldy et al ., 1985). Many of the processes 
involved in this interaction are also found in con- 
junction with auroral precipitation of energetic elec- 
trons (Gurnett et al .. 1986; Frank et al ., 1989). 

In the past, these processes have been studied in the 
laboratory as well as in space. The advantage of space 
experiments is the absence of walls and their con- 


taminating effects, which allows the study ot largcr- 
scalc phenomena such as the electromagnetic radi- 
ation from electron beams which usually have very 
long wavelengths compared with other scales in a 
plasma (Debye length. Larmor radius, etc ). On the 
other hand, a disadvantage is that it is more difficult 
to get detailed measurements. In the laboratory, the 
experiments can be repeated many times and measure- 
ments can be performed at different locations giving 
good temporal and spatial coverage. In space, this 
problem has been approached with the use of two or 
more separate payloads, one carrying a source and all 
payloads carrying a selection of plasma diagnostic 
instrumentation. Such configuration has proved very 
useful as demonstrated by the recent ECHO- 6 and -7 
sounding rocket experiments (Winckler et aL 1989; 
Erickson and Winckler. 1990) and the Spacelab-2 
flight of the space shuttle which released a plasma 
diagnostics package to free-fly during electron beam 
emission from the space shuttle orbiter (Gurnett et 
al ., 1986; Reeves et al., 1988a). 

The plasma disturbances are a consequence of the 
sudden ejection of large amounts of charge which 
create charge imbalances and large transient electric 
fields at beam turn on and off. In addition, the drift 
relative to the ambient plasma of the energetic primary 
beam electrons and the lower energy return current 
electrons and secondary electrons generate plasma 
waves in a broad frequency range. The hot plasma 
region around the beam flux-tube was investigated 
in the ECHO-1 experiment. Here three sub-payloads 
were injected upwards relative to the main payload at 
various angles to the magnetic field. The main payload 
carried an electron beam accelerator ejecting beams 
at a range of energies and currents (up to 36 keV, 
180 mA) as well as beam pitch angles. The turbulent 
region as observed from one of the sub-payloads, the 
Plasma Diagnostics Package ( PDP ), is shown in Fig. 
10 (Winckler et aL 1989). In the figure, the floating 
potential of the PDP is shown as function of per- 
pendicular distance to the beam flux-tube for a range 
of beam pitch angles (the floating potential is pro- 
portional to the temperature of the plasma). The 
dimension of the disturbed region generally increases 
with beam pitch angle. It is largest for 75 -85 pitch 
angle even though the beam in this case ejected down- 
wards and away from the PDP. This indicates that 
much of the disturbance is accounted for by secondary 
and scattered electrons in the flux-tube rather than 
the primary beam electrons. 

In the disturbed region around the beam flux-tube, 
low frequency electrostatic waves are excited in the 
range from essentially zero frequency and up to har- 
monics of the electron gyrofrequency (Shawhan et al ., 
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Fig 10 ECHO-1 observations of the hot plasma for three injection pitch angles in terms of the 

FLOATING POTENTIAL DECREASE OF THE PDP PAYLOAD BODY. PROPORTIONAL TO k T. AS A FUNCTION OF 
DISTANCE TRANSVERSE FROM THE BEAM INJECTION FIELD LINE. ALL FOR 36-kV, 180 mA INJECTIONS (FROM 

Winckler et al ., 1989). 


1984; Neubert et al 1986, 1991 ; Abe et ai, 1988; 
Frank et al ., 1989). The low frequency modes have 
been studied in the laboratory (Gekelman and Stenzel, 
1978; Stenzel, 1978a,b) and in computer simulations 
(Okuda and Ashour-Abdalla, 1988). The concensus 
is that the broad-banded noise below the lower hybrid 
frequency is ion acoustic noise generated by lower 
energy electrons flowing primarily as return currents 
to the payload. The disturbed region is, however, not 
symmetric around the magnetic field. Due to the 
motion of a spacecraft perpendicular to the field, there 
exists a turbulent wake region behind the beam 
(Frank et al., 1989; Reeves et al., 1990a). In the beam 
region, as in the wake, elevated electron temperatures 
and low frequency wave turbulence are found. The 
low frequency waves are then generated not just by 
return current electrons but, more generally, by lower 
energy electrons. 

Electromagnetic broad-banded noise in the whistler 
mode extending from the lower hybrid frequency and 
to the electron gyrofrequency is thought to be gen- 
erated by the primary beam electrons (Gurnett et al., 
1986). We return to this point in the following section. 

Localized primarily to the turbulent beam region 
are waves around the electron gyrofrequency and its 
harmonics, the plasma frequency and the upper 
hybrid frequency (Cartwright and Kellogg, 1974; 
Gurnett et al 1986; Neubert et al. , 1991). Computer 
simulations have shown that these high frequency 
wave modes are generated by the primary' beam elec- 
tron motion as well as by the induced azimuthal £ x B 
drift of beam electrons and ambient plasma electrons 


around the beam. The electric field here is either the 
radial field caused by the excess negative charge in 
the beam column or the electric field from waves in 
resonance with the particles (Winglee and Kellogg, 
1990 ; Okuda and Ashour-Abdalla, 1990). As a result, 
the beam electrons experience energy loss and scat- 
tering. 

Computer particle simulations have proven to be 
very useful in studying spacecraft charging, return 
current collection, and beam coherence. Using 2-D 
simulations (three velocity components) Winglee and 
Pritchett (1988) studied many aspects of the dynamics 
of electron beams emitted from spacecraft. In Fig. 1 1 
is shown the response of the plasma to the ejection of 
an electron beam at a time at pb t = 120, w here <o pb / is 
the electron beam plasma frequency. The magnetic 
field is along the x-axis and the beam is emitted at a 
finite angle to B. Figure 1 la shows contours of the 
charge density of the electron beam and the spacecraft, 
which is modelled as a planar conductor. The helical 
primary beam can be seen for about two gyrations, 
embedded in beam electrons scattered by electrostatic 
fields from the beam space charge, plasma wave noise, 
and the spacecraft electric fields. Figure 1 lb shows the 
ambient plasma electron density. The shaded regions 
are those with densities smaller than the unperturbed 
density. The depletion is brought about by the col- 
lection of return currents to the spacecraft as shown 
by the arrows. Behind the spacecraft this region is of 
a dimension perpendicular to B which is comparable 
with the spacecraft dimension. In front of the space- 
craft, the region is surrounding the primary beam 
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Fig 11 Particle simulation of the injection of electron beams from spacecraft. 

Contours of (a) the charge density of the beam and the spacecraft, (b) the plasma electron density, and 
(c) the plasma ion density at = 120. Successive contours in (a) differ by a factor of 3.16 while those 
in (b) and (c) decrease linearly from the maximum value down to 8% of the maximum value. The arrows 
in (b) and (c) indicate the direction of the plasma flows in response to the beam injection. The shaded 
areas indicate the region where the density is smaller than in the initial plasma density (from Winglee and 

Pritchett, 1988). 
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column. Current closure is in this simulation achieved 
by the plasma ion motion shown in Fig. lie. which 
allows ions to flow across the magnetic field to a 
distance limited by their Larmor radius. 

The beam structure obtained in the 2-D simulation, 
is consistent with optical observations of the electron 
beam structure during the STS- 3 mission ol the space 
shuttle (Banks and Raitt, 1988). It was found that 
the electron beam initially moved outwards along a 
classical helical trajectory. However, within a short 
distance electrons were deflected into paths parallel 
and anti-parallel to the magnetic field lines threading 
the primary helix. This resulted in the formation of a 
thin-walled, hot. electron torus with a radius equal to 
the primary beam gyroradius. 

Winglee and Pritchett (1988) also found that a 
build-up of negative space-charge in the beam region 
close to the spacecraft may occur if the beam stag- 
nation time f s is smaller than the plasma response time 
/ rp . For a given beam current, the ratio of the two 
time-scales depends on the beam density, the beam 
energy and the area of the spacecraft. The negative 
space-charge region is often referred to as a virtual 
cathode. The formation of a virtual cathode partially 
inhibits beam escape, destroys the spatial coherence 
of the beam, and creates additional plasma wave tur- 
bulence (Pritchett, 1990). 

Early 1- and 2-D simulations gave rather restrictive 
values for the beam stagnation time and questions 
arose about beam escape from spacecraft in iono- 
spheric beam experiments (Pntchett and Winglee, 
1987; Winglee and Pritchett, 1988). However, other 
studies argued that experimental observations indi- 
cated beam escape, or at least that a substantial por- 
tion of the beam escaped the spacecraft (Neubert et 
tv/., 1988 ; Farrell, 1990). The discrepancy was resolved 
when 3-D simulations were performed. In three 
dimensions the criteria for beam escape are less restric- 
tive because electrostatic fields decrease more rapidly 
with distance than in two dimensions. Simulations 
and experimental observations are now in agreement 
on this point : namely, that in most cases in past exper- 
iments, the beams were able to leave the payload 
environment (Okuda and Ashour-Abdallah, 1991 ; 
Pntchett, 1991). 

Nevertheless, the virtual cathode effect is still 
important for electron beam emissions in space, and 
while it may not have limited the beam escape in past 
experiments, it certainly affects beam coherence. In 
addition, it is probably responsible for much of the 
electrostatic noise observed from the beam-emitting 
spacecraft both at low frequencies as observed in the 
SEPAC space shuttle experiment (Neubert et a/., 
1986; Cai et al. . 1987) and at high frequencies as 


suggested by Neubert et al. (1991). Furthermore, 
return current electrons with energies in excess of the 
beam accelerator energy observed in the past (Man- 
agadze et al ., 1988; Wiihelm et aL 1984; Waterman 
et al .. 1988) and optical emissions (Sasaki et al 
1985a) are probably a result of such space-charge 
effects in the near-environment of the spacecraft 
rather than caused by a beam-plasma instability or a 
beam-plasma discharge. 

Transient excursions of the spacecraft potential to 
levels well above the beam energy, often called “super- 
charging". have occasionally been observed (Sagdeev 
et al ., 1981; Maehlum et aL 1988). It has been 
suggested that these high fields are associated with the 
formation of a virtual cathode. It is proposed that 
electrons are accelerated to high energies in the cath- 
ode region and that a percentage of these escape at 
energies which exceed the beam energy thereby driving 
the spacecraft to large positive potentials (Maehlum 
et al 1988; Managadze et aL 1983, 1988; Denig et 
aL 1990). 

The results on plasma heating, wave excitation, and 
beam dynamics can be summarized as : 

(7) The 3-D structure of the disturbed electron 
beam region is now' explored, revealing the dimension 
of the disturbed region to be about 4 beam electron 
Larmor radii perpendicular to the magnetic field. It 
has also been established that there exists a turbulent 
wake region behind the beam flux-tube. 

(8) It is clear that transients in electrostatic poten- 
tials as well as return currents, azimuthal currents 
around the beam, and the primary beam current are 
important for the generation of broad-band radiation 
observed during electron beam ejection. 

(9) Beam dynamics including electrostatic beam 
expansion, return current structure, and the phenom- 
enon of virtual cathode and other space charge effects 
are now better understood. 

(10) The success of the past decade in the under- 
standing of the complete beam spacecraft system is 
due to the synergy of detailed observ ations from mul- 
tiple platforms, analytical models, and computer par- 
ticle simulations. 

2.4. Neutral gas emissions 

The normal electrical neutralization of spacecraft 
at high altitudes relies solely on the ability to collect 
sufficient current from the ambient plasma since ion- 
ization of the neutral atmosphere provides little 
additional plasma. While the beam-system inter- 
actions mentioned previously have been thought to 
increase the electron current across the magnetic field, 
and thereby increase the return current through wave- 
induced particle diffusion (Linson. 1969; Hastings, 
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198^). there are practical limits to the efficiency of such 
processes. Furthermore, while laboratory experiments 
indicate that such processes occur in this environment, 
it has not been possible to establish their importance 
in space experiments. 

Another possibility for increasing the return current 
to a beam-emitting payload is through ionization of 
pavload gases introduced accidentally through out- 
gassing or deliberately, as in the case of attitude con- 
trol thrusters. For example, in the ECHO - 4 rocket 
experiment it was found that the injection of N 2 gas 
by attitude control thrusters on the beam-emitting 
payload had the effect of lowering the spacecraft poten- 
tial during beam emissions (Israelson and Winckler, 
1979). This was also observed in the SEP AC 
experiment (Sasaki et al . , 1985b) as well as in the 
CHARGE-2 experiment. In the case of CHARGE-2, 
N : emissions from the beam-emitting mother as 
well as from the tethered, but passive daughter, 
reduced the spacecraft potential (Gilchrist et al ., 1990 ; 
Banks et aL 1990). The reduction in potential is 
thought to be brought about by the ionization of the 
neutral gas cloud. The ionization is not dependent on 
the primary beam electrons but is triggered by return 
current electrons energized in the sheath. This con- 
clusion was reached from the fact that the CHARGE- 
2 daughter payload was located so far from the mother 
payload that the gas cloud could not have reached the 
beam region in the time it took the daughter to drop 
in potential. As a result of the decreased potential, 
wave turbulence was decreased (Neubert et al.. 1991 ), 
and the beam escaped with a larger fraction of its 
energy as inferred from ground-based radar obser- 
vations of the plasma density enhancement caused by 
BAI (Gilchrist. 1991). The reduction of spacecraft 
potential and easier escape of the beam from the pay- 
load in the presence of ambient or thruster neutral 
gases have also been found in computer simulations 
(Winglee, 1990). 

Another means of controlling the spacecraft poten- 
tial involves the use of the so-called plasma contactor, 
or hollow cathode. It operates by creating a dense 
plasma cloud in the vicinity of the spacecraft which 
then acts as a bridge between the more dilute back- 
ground plasma and the spacecraft (Patterson. 1987; 
Williams et al.. 1987; Katz and Davis, 1987). The 
plasma contactor in its normal mode of operation 
relies on the ionization within the hollow cathode of 
neutral gas flowing out from a pressure chamber. At 
positive bias potential, the contactor may operate in 
the so-called ignited mode in which the ionization 
occurs in the gas plume in the sheath region. Thus, 
the ignited mode is simply the ejection of neutral gas 
as observed in the CHARGE - 2 experiment. The 


importance of this discovery is that with this method 
it is fairly simple and inexpensive to control large 
positive potentials and that such control seems to 
reduce plasma wave noise and allow the beam to 
escape with almost its complete accelerator energy. 
The limits in terms of beam current that can be neu- 
tralized in this way and the gas densities and flow rates 
required need to be explored in future experiments. 

The results given above can be summarized : 

(11) Large, positive payload potentials can be 
reduced by the ejection and subsequent ionization of 
neutral gas. The limits of beam currents that can be 
ejected during neutral gas emissions and the gas den- 
sity and flow rates required to neutralize the payload 
need to be quantified with accurate experiments in 
space. 

(12) When the payload potential is reduced by neu- 
tral gas emissions, the plasma wave noise is reduced 
and the beam escapes with a larger fraction of its 
accelerator energy. 

3. WHISTLER MODE RADIATION FROM ELECTRON 
BEAMS 

3.F Background 

We now explore some basic features of wave gen- 
eration from pulsed and d.c. electron beam exper- 
iments in space plasmas. We pay particular attention 
to whistler wave generation in ionospheric plasmas 
and to the limits on the electron beam energy for wave 
generation in the two principal resonances, Landau 
resonance and cyclotron resonance. The standard 
way of calculating the radiation from a single particle 
spiralling in a magnetized cold background plasma is 
to solve Maxwell's equations with the moving particle 
as the current source (Mansfield, 1967; McKenzie, 
1967). In this treatment, the radiation is found to be 
particularly intense for waves satisfying the resonance 
condition ; 

R = — tu = 0. (13) 

In equation (13), ca* is the relativistic beam particle 
cyclotron frequency : 

ay* = tu c (l — v 2 /c 2 ) 1 ", (14) 

where r is the particle velocity and c is the velocity of 
light. The component of the wave vector parallel to 
the magnetic field, k ,, is given by 

k t = k cos 0, (15) 

where k is the wavenumber and 0 the angle of the 
wavenumber to the magnetic field. Similarly, for the 
particle velocity along the magnetic field : 
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r = r cos a. (16) 

where y. is the pitch angle ot the particle. The wave 
frequency o) is related to k through the plasma dis- 
persion relation. 

The parameter .5 takes on the value 0 for Landau 
resonance where Cherenkov radiation is emitted, 1 
for cyclotron resonance, and — 1 for anomalous cyclo- 
tron resonance. Larger integer values denote higher 
order cyclotron resonances. Condition (13) is often 
referred to as the first-order resonance. Since the back- 
ground plasma is usually inhomogeneous, waves and 
particles in resonance at one point will become off- 
resonant as they propagate away to other regions. 
Enhanced interaction is achieved when waves and 
particles are also in second-order resonance, namely 
w hen : 

dR( r)/dr = 0, (17) 

where z is the coordinate along the background mag- 
netic field (Neubert el al.. 1987). 

As an example, we consider the case of whistler 
mode radiation. An approximate expression for the 
refractive index g for to > w LHR and oj « < 10 ^ is : 

fi = CDpe [a>(to cc cos 0 — oj)] ’ ", (18) 

where g — kc oj . is the electron plasma frequency, 
and u> ce the electron gyrofrequency. The resonance 
cone angle is the angle at which the refractive index 
and the wavenumber goes to infinity : 

cos 0 T = (19) 

Real solutions to equation (18) require 0 < 0 T . Figure 
12 shows a sketch of the whistler mode refractive 
index surface as function of 0 for two values of the 
wave frequency. Also shown is the direction of the 
group velocity, which is perpendicular to the surface. 



Fig. 12. The whistler mode refractive index surface. 
The arrows indicate the direction of the group velocity, which 
is perpendicular to the surface. 


Foro < o), c 2 (satisfied by w, in Fig. 12) the refractive 
index surface has a dimple around 0 = 0 such that a 
minimum in g occurs for 0 - 0 t . the so-called Gen- 
drin angle. For oj > the dimple disappears (as 
is the case for oj z in Fig. 12) and the minimum in ji, 
occurs for 0 = 0. The minimum value of ji, found 
from equation (18) is : 

/W = 2to pt /w a: («<w„) (20) 

and the Gendrin angle is : 

cos = lojiajct {oj < ctfce/2). (21) 

The relativistic expression for the electron energy. 

normalized to the electron rest mass (511 keV) is : 

£ = (1 —v 2 jc 2 ) - 1 2 — 1, (22) 

where the resonance beam electron velocity found 
from equation (13) is: 

vjc = (23) 

with g* = l — jcuJ oj. 

From equation (23) it is seen that at a given wave 
frequency there is a maximum beam energy for which 
resonance is possible. This energy corresponds to the 
minimum in / 1 ,. A minimum of the beam energy can 
be found by noting that small wavelength whistler 
waves have a very large refractive index and are pri- 
marily electrostatic in nature. As a result they are 
strongly damped as they propagate through a plasma. 
Using a wavelength of I m as a practical lower limit 
for whistler mode propagation, this corresponds to a 
maximum of k = 2 tt m " '. 

The maximum and minimum electron energies for 
Landau resonance are shown in Fig. 13 and for cyclo- 
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Fig 13. Maximum and minimum resonance energy for 
Cherenkov resonance 
Shaded areas indicate regions of resonance. 
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iron resonance in Fig. 14. The region allowing res- 
onance with w holler mode waves is marked voih I 
(Region II is discussed later.) Also shown is a curve 
through this region which corresponds to the res- 
onance energy for 0 — () . It is assumed that o 2tt 
and rj,,, 2 Jr are 10'' s 1 and 5 .< 1 ()*’ s respectively, 
which are values representative tor the peak of the 
ionosphere. The ion population is exclusively O * . The 
figures show the resonance energies according to the 
Sux ( 1962) formulation of the cold plasma dispersion 
relation and includes relativistic effects on the beam 
electrons. It is assumed that the pilch angle, y. = 0 . 

At frequencies below one half the electron gyro- 
frequency and above the lower hybrid frequency. 
/( Mk . the maximum energy for Landau resonance is 
approximately constant in wave frequency for given 
background plasma parameters [from equations (20) 
and (23)]. In the ionosphere the ambient plasma den- 
sity and therefore the plasma frequencies may have 
large fluctuations depending on local time, geographic 
location, and even small variations in altitude. As a 
result, the maximum energy for Landau resonance 
may also vary and will generally increase with decreas- 
ing density. Many electron beam experiments per- 
formed in the past have been in the range from 1 to 5 
keV and were therefore expected to generate waves 
through Landau resonance in the whistler mode. Some 
experiments, in particular the ECHO scries, used 
higher energies, up to 40 keV. and the degree to which 
Cherenkov radiation was generated via Landau res- 
onance depended on the spacecraft altitude and iono- 
spheric conditions. 


Cyclotron Resonance 



Frequency (Hz) 

Liti. 14. Maximim a no minimi m risonaml iniroy ior 

CYCLOTRON Rl-SONANd 
Shaded areas indicate regions of resonance. 


Lor low frequency waves in cyclotron resonance. </, 
becomes large and. consequently, the resonance 
energy becomes large, as can be seen Irom equation 
(23). As the wave frequency approaches the lower 
hybrid frequency, the resonance cone angle 
approaches 90 and the maximum as well as the mini- 
mum resonance particle energy go to infinity. Because 
of the higher energy required for cyclotron resonance, 
this resonance has not been studied to the same degree 
as Cherenkov resonance. Future experiments with 
higher energy accelerators may study cyclotron res- 
onance. It will he relatively simple to distinguish 
between cyclotron and Landau resonance. In the 
latter, the wave energy propagates in the same direc- 
tion as the beam electrons, while in the former the 
waves propagate opposite to the beam electrons. 

Below the lower hybrid resonance the topology of 
the wave normal surface changes. The resonance cone 
disappears and the wave normal surface becomes 
closed. The maximum in energy found for frequencies 
above the lower hybrid resonance becomes a mini- 
mum in energy below r the lower hybrid resonance. The 
region of electron energies in resonance here is marked 
II in Figs 1 3 and 14. 

3.2. Radiation from J.e. electron beams 

On auroral field lines charged particles are accel- 
erated by localized electric fields thought to be located 
at altitudes of 2-4 Earth radii. The result is the cre- 
ation of an electron plasma with beam-like dis- 
tributions in velocity. These electrons excite aurora 
as they penetrate the upper atmosphere. They also 
generate plasma waves over a wide frequency range. 
Figure 15 shows a typical frequency spectrum of the 
electric field intensity observed by the DE 1 satellite as 
it traverses the auroral oval from the polar cap and 
to the plasmasphere (Gurnett et a /.. 19X3). Above the 
electron gy rofrequcncy auroral kilometric radiation 
(AKR) and emissions at the upper hybrid frequency 
are seen, while below the gy rofrequcncy. in the whis- 
tler mode frequency range, auroral hiss and plasma- 
spheric hiss are seen as the spacecraft enters the 
plasmasphere. In particular we point out the 
funnel shape of the auroral hiss as it appears on a 
frequency-time plot as shown in Fig 15. This shape 
is consistent w ith the idea that the whistler wave noise 
is generated with wave normals close to the resonance 
cone through Landau resonance with energetic aur- 
oral electrons accelerated at some localized altitude 
below the spacecraft. 

The funnel shape is a result of the propagation 
characteristics of whistler mode waves at the res- 
onance cone. The higher the wave frequency, the nar- 
rower the resonance cone; because the direction of 
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the group \clociu lor wa\c^ at the resonanee cone F 
perpendicuiar to k. the angle of propagation to H 
increase^ with wave lrequenc> . Thus high-frequenev 
emission-* are first encountered as the spacecraft 
approaches the field line. In the case of a point source. 
onl> a V-shaped outline of the lunnel will be seen. 
Such emissions are frequently observed and have been 
named "saucers” (James. 1976). A "tilled-in tunnel 
requires a beam (or sheath) source which extends 
from a given point in space. 

Emissions of a similar nature have been observed 
in artificial beam experiments. In the Spaeelah-2 exper- 
iments made with the space shuttle, a small satellite, 
the Plasma Diagnostics Package {PDP). was released 
from the shuttle and separated to distances of about 
300 m during a 6 h period. During that time, the shut- 
tle was maneuvered to provide four magnetic con- 
junctions of the PDP and the shuttle. During two of 
these, an electron beam (with electron energies of 1 
keV. and currents of 50. 100. or 150 mA) was emitted 
from the shuttle. The electric field intensity observed 
by the PDP during the ejection of a continuous beam 
at 50 mA is shown in Fig. 16 (Gurnett vt al. . 1986). 
Over frequencies between about 30 kHz and up to the 
electron cyclotron frequency a funnel-shaped emis- 
sion is seen which bears close resemblance to the 
auroral hiss shown in Fig. 15. As the PDP ap- 
proaches the beam flux-tube the high frequency emis- 
sions are observed first followed by lower frequency 
emissions. When the PDP crosses the beam, emissions 
are also observed close to the local upper hybrid fre- 
quency inside the beam. 

The similarities of the natural and beam-generated 
emission led to the conclusion that they were a result 
of the same process. In the case of the Spacvlah- 2 
beam, it was found that the spectral shape on lime 
frequency plots as well as the electric field polarization 
were consistent with the assumption of Landau res- 
onance between beam electrons and whistler waves at 
the resonance cone (Farrell et al .. 1988). 

The Spacvlah- 2 observations allowed the quan- 
titative comparison of theories for wave generation 
by electron beams with the observations. Studies in 
the past of natural!) occurring plasma waves and par- 
ticle distributions using the formulation of Mansfield 
(1967) and McKenzie (1967) led to the conclusion 
that incoherent Cherenkov radiation (the radiation 
from each electron is added incoherently) was several 
orders of magnitude lower than observed (Taylor and 
Shaw han. 1974). Because the power radiated by 
coherent radiation is proportional to .V\ where V is 
the number of particles, while the power from inco- 
herent radiation is proportional to A’, an assumption 
of some coherency will increase the power and bring 


it in agreement with observations. Similar conclusions 
were reached for the Spacclah-2 electron beam (Farrell 
et al.. 1988). It was found that the radiation from the 
beam was orders of magnitude above expectations tor 
an incoherent process and that some coherency must 
have been present. Following an idea of Bell ( 1968). 
analytical calculations and l-D computer simulations 
were presented by Farrell et al. (1989). which dem- 
onstrated that density perturbations in the electron 
beam density, caused by Langmuir oscillations in the 
beam, as observed during magnetic conjunction, can 
provide the necessary coherence. The density per- 
turbations were found to propagate with a velocity- 
close to the beam electron velocity and could therefore 
radiate through the Landau resonance. 

Later simulation studies in two dimensions con- 
firmed this conclusion with minor modifications. Prit- 
chett et al (1989) found that the beam-plasma inter- 
action creates a current structure that radiates much 
like an antenna through the Cherenkov radiation 
mechanism. In their study, the beam region carried 
no net current at the start of the simulation. This 
was achieved by giving the higher density background 
plasma an appropriate low energy drift velocity in the 
opposite direction to the beam velocity. Although the 
average net current is zero, radiation is still generated 
because of the electron current perturbations created 
in the beam bv micro-instabilities. 

Other studies of electron beam radiation include 
that of Wong and Lin (1990) which presents an ana- 
lytical model for the modes generated inside the beam 
and the coupling at the beam edge of these modes to 
the ambient plasma. It was found lhal this process 
could generate whistler mode radiation which propa- 
gates in the ambient plasma with k close to the res- 
onance cone as observed in the Spacvlah- 2 exper- 
iments. A particle simulation of the Spaeelah-2 
experimental results was presented in Omura and 
Matsumoto (1988). The simulation was done in two 
dimensions with periodic boundary conditions in the 
dimension along the beam. Such boundary conditions 
amount to simulating an infinite beam. Even so. the 
results showed wave spectra as a function of distance 
from the beam which resembled the funnel-shaped 
spectra oi' Spaeelah-2. They suggested that the funnel 
was a result of frequency-dependent damping as the 
waves propagate away from the beam. Furthermore, 
they found that the wavelength of the electrostatic 
noise Generated inside the beam column did not match 
the wavelength of the radiated whistler wave noise 
and concluded that the radiation was not caused by 
bunching of the beam electrons. 

The results of radiation from continuous electron 
beams can be summarized as follows : 
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ELECTRIC FIELD, DE-I, OCTOBER 15, DAY 288, 1981 

Fig ISA representative spectrogram of the electric field INTENSITIES FOR a sight-side crossing 
OF THE AURORAL FIELD LINES AS OBSERVED BY DE- 1 (FROM GlRNETT el al.. 1983). 


171 



PDP, DAY 213, AUG. I, 1985 

Fl(i. 16. A SPH TRIXiRAM OF THE ELECTRIC HI LO INTENSITY DURING ELECTRON BEAM FJECTIONS FROM THF 

SPACE SHUTTLE 

The free-flying PDP is in magnetic conjunction with the shuttle and in the beam region at the time indicated 

by the arrow (from Gurnett el al., 19X6). 
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Fig. 17. The magnetic field in the frequency ranges 0 10. 20 10. AND 20 30 kHz AS a function of 

time:. 

The PDP is free-flying and the FPEG is pulsed at 1 .22 kHz (from Bush et aL 1987; Neubert cl ai. 1988). 
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(13) The consensus from experiments a nil eompuier 
simulations is that in past experiments, employing 
kiloelectron\o!t energy electron beams, whistler wave 
radiation is generated through Cherenkov resonance 
with beam electrons. 

(14) The intensity of the radiation is between the 
levels expected for incoherent and coherent radiation. 
Micro-instabilities in the beam are found to be impor- 
tant for the establishment of quasi-coherency and 
thereby for the creation of the structure of the radi- 
ating current. 

(15) The details of the formation of the radiating 
current structure are not well understood. It is thought 
that Langmuir oscillations play an important role in 
bunching the electrons. Spacecraft potential oscil- 
lations may also be important. 

3.3. Radiation from pulsed electron beams 

A pulsed electron beam generates electromagnetic 
radiation at the pulsing frequency and its harmonics. 
The radiation signature is therefore easier to dis- 
tinguish from ambient plasma noise than the radiation 
from continuous beams. Furthermore, the spectrum 
of the radiation contains information on the spatial 
structure of the pulsed beam which can be compared 
with the format of the injected beam. In this way 
additional detailed information on beam dynamics 
and radiation can be obtained. Since the very first 
space-borne electron beam experiment (Hess et a /.. 
1971), it has been conjectured that pulsed electron 
beams could be used as virtual antennas lor com- 
munication purposes, for instance between space and 
ground or between two platforms in space. Beams 
have an advantage over conventional electric dipole 
antennas or magnetic dipole loop antennas, which 
is particularly important in the VLF whistler mode 
frequency range; because ol the large wavelength of 
VLF waves it takes very long electric dipoles or large 
loops to efficiently radiate VLF waves. 

The STS-} mission of the space shuttle orbiler, 
flown in 1982. carried an electron beam accelerator 
especially designed to study the radiation from pulsed 
electron beams. The accelerator, the Fast Pulse Elec- 
tron Generator (FPEG) emitted l keV energy electron 
beams with currents ol 50, 100 or 150 mA. The beam 
could be square-wave modulated from d.c. and up to 
several hundred kilohertz w ith a pulse rise-time of the 
order of 10 *s (Banks et al .. 1987). The elec- 
tromagnetic radiation generated by beams pulsed at 
ELF and VLF frequencies was observed with wave 
receivers mounted on a PDP , w hich was either stored 
in the payload bay. or mounted on the Remote Man- 
ipulating System (RMS) (Shawhun et al.. 1984; 
Reeves et al .. 1988b). 


Attempts were made to detect the radiation at 
remote locations m the magnetosphere by executing 
beam emission sequences coordinated with the Dt I 
satellite at times when the shuttle and the satellite were 
in magnetic conjunction. The results were negative; 
however it was later shown that either the STS-} 
altitude was such that the electrons hit the main body 
of the orbiter. or ra> tracings showed that it was not 
possible for whistler mode waves to propagate from 
the STS- 3 location up to the vicinity of the DE 1 
satellite (Inan etal.. 1984). The FPEG PDP was flown 
again in the Spacelah- 2 mission of 1985 as mentioned 
earlier. Radiation was observed on the PDP free- 
flycr out to its maximum separation distance of about 
300 m (Bush et al .. 1 987 ; Reeves et al.. 1 988a ; Neubert 
et al ., 1988). Again no radiation was observed on the 
ground or on DE 1 . 

Results from other experiments studying the radi- 
ation characteristics of pulsed electron beam emis- 
sions are reported in Cartwright and Kellogg (1974) ; 
Reme (1980); Holzworth and Koons (1981); 
Winckler et al. (1984. 1985); Kellogg et al. (1986); 
Sasaki et al. (1987, 1988); Gocrke et al. (1990); 
Neubert et al. (1991). In most of these experiments, 
attempts were made to detect the wave fields on the 
ground. The results to date have proved negative, and 
the largest distance from the beam source that VLF 
wave fields have been observed is of the order ol 1-2 
km. 

In parallel with the experimental efforts involving 
the FPEG. a theory for the radiation from the square- 
wave modulated electron beam was developed. The 
first model considered the coherent radiation Irom a 
finite train of pulses in jected into a magnetized plasma 
(Harkerand Banks. 1983). For simplicity, each pulse 
was assumed to traverse an unbounded path length 
( — x to + /.), The beam was assumed to radiate 
coherently and the radiation was determined for the 
frequency range from the lower hybrid frequency to 
the electron gyrofrcquency. In Harker and Banks 
(1985) the power radiated by a beam propagating into 
a half space (0 to 4- x ) w as calculated for frequencies 
both below r and above the lower hybrid frequency. 
The power was found for the tar-field region, which 
is the region at distances beyond a few perpendicular 
wavelengths from the beam source. Since the Space- 
lah- 2 observations were performed in the near-field 
region a near-field model was developed (Harker and 
Banks. 1987) giving the complete polarization of the 
electric field. Observations showed, however, that 
associated with the pulsed beam injection, con- 
siderable broad-band electrostatic noise was gen- 
erated just as for the case of continuous beam injec- 
tion. This made it difficult at times to determine the 
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spectral component on the pulsing frequency or its 
harmonics. The magnetic field, on the other hand, was 
very clear and undisturbed from broad-hand noise. 
The calculations were therefore extended to include 
the near-Held components of the magnetic held (Ncu- 
bert and Marker, 1988). Finally. the theory has been 
extended to allow calculations of signal levels on the 
ground from pulsed beam emissions in space (Harker 
et al.. 1991). 

The theory of Harker and Banks is based on the 
work of McKenzie (1967). It assumes that the radi- 
ating current is that ol an infinitely thin electron beam 
spiralling around the geomagnetic field under the 
influence of the Lorentz force. The beam is square- 
wave modulated (pulsed) and the beam electrons 
within a pulse arc assumed to radiate coherently. 
Effects of electrostatic expansion of the helical beam, 
return currents and other processes which destroy the 
ideal beam structure are accounted for through an 
attenuation factor which assumes that the ability of 
the beam to radiate coherently decreases expo- 
nentially with distance from the beam accelerator. 

The Harker and Banks theory was developed par- 
ticularly for the square-wave modulated electron 
beam of the F PEG. Other treatments include the radi- 
ation from pulsed electron beams at zero degree pitch 
angle (Lavergnat and Lehner. 1984; Ohnuki and 
Aduchi. 1984) and sinusoidally modulated electron 
beams at arbitrary pitch angles (Lavergnat et a /.. 
1984). 

As a result of the assumption of a square-wave 
modulated current source, the radiated fields contain 
a factor D. the so-called duty cycle factor; 

D = 1 r/ sin d). (24) 

where bd is the duty cycle, b is the on-time for a pulse, 
and d is the total modulation period. The parameter rj 
is the harmonic number. The factor D arises naturally 
from the Fourier transform of a square-wave tunc- 
tion. 

During the Spaeelab- 2 experiment, several pre- 
dictions of the Harker and Banks theory were tested. 
The PDP t ree-flight provided observations of the elec- 
tric and magnetic fields during a beam sequence in 
which the beam (100 mA) was pulsed with 50% duty 
cycle at 1.22 kHz. The PDP was located during this 
beam sequence at distances from a few meters to 1 50 m 
from the electron beam. A spectrogram of the mag- 
netic field is shown in Fig. 17 (Bush et a /., 1987; 
Neubert et al ., 1988). The relative signal intensity is 
color coded showing the fundamental and the odd 
harmonics of the modulation frequency as horizontal 
lines. In the first half of the panel (0-25 s) the fre- 
quency band is from 0 to 1 0 kHz. the following quarter 


(25 38 s). the band is 20 10 kHz with 20 kHz at the 
bottom of the frequency scale, and in the last quarter 
the band is from 20 to 30 kHz. During the time interval 
shown the distance of the PDP was 75 m per- 
pendicular from the beam. 

The Fourier transform of a square-wave function 
with a 50° o duly cycle has vanishing even harmonics 
and odd harmonics with amplitudes that varies as 1 f. 
Figure 18 shows the magnetic field amplitude of the 
narrow-band emissions observed at three different 
locations of the PDP relative to the beam (Reeves et 
al ., 1988a). As can be seen, the amplitude does 
decrease approximately as 1 /. Note also that even 
harmonics arc present, although weaker by an order 
of magnitude. This indicates that the beam is losing 
some of its square-wave form within the 300m dis- 
tance of the PDP and the shuttle. However, as pointed 
out by Reeves et al. (1990b), the duly cycle factor is 
very sensitive to changes in the duty cycle. A change 
from 50 to 50.5% is enough to account for the ampli- 
tude of the even harmonics. 

The beam resonance condition determines k v The 
dispersion relation has two roots or two values of k 
for given values of k and to. In certain frequency 
bands one of these roots becomes imaginary, which 
corresponds to evanescent waves. In the near-field 
these waves can give important contributions to the 
electromagnetic field. In the theory ol Harker and 
Banks root! is evanescent below / lHR while root2 is 
continuous across / LH r Above /j HR root l cor- 
responds to waves generated at the resonance cone, 
while root2 corresponds to waves generated in the 
central hump of the refractive index surface. Studies 
of the amplitude of the harmonics as a function of 
distance trom the beam show that the signal amplitude 
roughly follows the predictions for the root2 waves. 
For the Spacclah-2 frce-flight./ LHR was around 3 kHz. 
and therefore the first harmonic was below / LHR while 
the third and higher harmonics were above / LHR . The 
root 1 solutions for the higher harmonics have, accord- 
ing to the theory, much larger amplitudes than 
observed. Furthermore, the first harmonic (which is 
evanescent lor root! ) should be strong very close to 
the beam, and decrease rapidly with distance until the 
rool2 waves become dominant. This behavior was 
not observed in the magnetic field data. Rather, the 
amplitude varied with distance much as predicted for 
root2 waves (Reeves et al.. 1990a). 

As shown in Fig. 19, the first harmonic of the elec- 
tric field below / LH r did decrease relatively rapidly 
with distance ( Reeves et al.. 1988a). At the same time, 
emissions around the second and third harmonics, 
broad-band as well as narrow-band, became relatively 
dominant. It is possible that these waves are lower 
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Fig 18. Magnetic field amplitude of the harmonics of the beam pulsing frequency <1-22 kHz) at 

THREE DIFFERENT LOCATIONS OF THE PDP RELATIVE TO THE BEAM (FROM REEVES el ui, 1988a). 


hybrid waves as suggested by computer particle simu- 
lations (Matsumoio and Fukuchi, 1985; Hwang and 
Okuda, 1989). 

The dependence of the wave field amplitude on the 


duty cycle has been studied by injecting a sequence of 
pulses with constant pulsing frequency and varying 
the duty cycle. Figure 20 shows the magnetic field 
amplitude observed in the payload bay of the first 
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Fig. 19. Electric field amplitude of the harmonics of the beam pulsing frequency (1.22 kHz) at 
DIFFERENT LOCATIONS OF THE PDP (FROM REEVES el u/.. 1988a). 
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The solid line is the duly cycle factor D arbitrarily normalized 
to 0.1 (from Reeves ci al.. 1 990b i. 


harmonic as a function of duty cycle (Reeves a a!., 
1 990b The factor D is normalized to the arbitrary 
value of 0.1 in order to enable the comparison with the 
observations. As can be seen, the agreement between 
observations and predictions is quite good. As the 
duty cycle goes to zero, the average current goes to 
zero and the field amplitude approaches some small 
background level. For the duly cycle approaching 
100%, the beam becomes continuous and only broad* 
band noise is generated. Here, the amplitude again 
approaches a very small value. One question that has 
remained unanswered in the past has been whether 
the sudden onset and turn-off of the beam, that is. the 
leading or trailing edge of the beam or sudden changes 
in the spacecraft potential were responsible for the 
radiation. In this case, radiation should be inde- 
pendent of the duty cycle. Thus, the Spaccluh-2 study 
shows that the source at least to some approximation 
is square-wave modulated as predicted by theory . 

The polarization of the beam-generated fields was 
studied in Ncuberl ci al. (1990c). Although only one 
component of the electric and one component of the 
magnetic field was measured, the spin of the POP was 
utilized to find the polarization of the fields in the spin 
plane. It was found that the observed polarization did 
not agree well with the predicted polarization and 
that therefore the current structure departed from 
the helical structure assumed in the model. Another 
possibility comes from the fact that the theory pre- 
sented in Harker and Banks (1987) neglects the con- 
tribution from a branch cut in the integration over 
k which corresponds to ray propagation, and only 
considers the simple poles which correspond to sur- 
face waves. The model is, therefore, incomplete as it 
stands. 

A final comparison with theory is shown in Fig. 21 
(Reeves ct al.. 1990b). The amplitude of the fun- 
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Predicted curves are >hown for Cherenkov root I and root2. 
The size of the markers indicate the duty cycle (from Reeves 
el al.. 1990b). 


damental component is plotted as function of 
frequency. This coverage in frequency required a 
compromise on the duty cycle. The duty cycle used is 
shown by the size of the markers. 50% duty cycle 
having the largest marker. At this time the POP was 
located in the payload bay. Also shown arc the pre- 
dictions of the model for rootl and root2. Again it 
can be concluded that the model predictions for root- 
are in general agreement with observations while 
rootl predicts amplitudes many orders of magnitude 
above observations. 

The results from pulsed electron beam studies can 
be summarized as follows : 

(16) Observations show that the spectral structure 
of the wave emissions closely follows the spectral 
structure of the injected beam current. A square-wave 
modulated beam generates harmonics as expected 
from the Fourier transform of a square-wave func- 
tion. 

( 17) A closer comparison with the theory of Harker 
and Banks (198 7 ) reveals that the amplitude depen- 
dence of individual harmonic components with fre- 
quency or distance from the beam is not as expected 
from theory. In addition, the observed polarization is 
not as expected. The amplitude is in general agreement 
with the predictions of root2. 

( 18) A conclusion must be that the model does not 
assume the correct structure of the radiating current. 
For instance, a current structure as expected from 
simulations of continuous beam emissions, when 
modulated by a square wave, will give the observed 
spectral structure. The amplitude and polarization 
from such a model may agree belter with obser- 
vations. 
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4 DISCI SSIOS 

I jom the 1* Nummary points given throughout t hi > 
paper n is seen that a basic understanding of many 
processes associated with electron beam ejection in 
space plasmas has been achieved. This knowledge has 
emerged from a synthesis of experiments, analytical 
modelline, and computer simulations. As a result, 
beam ejection can now be controlled to a higher degree 
than before, allowing the ejection of larger currents, 
increasing the beam coherence and the escape energy 
of the beam, and reducing the plasma wave noise. 
These are crucial skills for some of the applications of 
electron beams to improv ing understanding of natural 
processes of plasma wave generation. 

Experiments in the future are needed to study still 
unresolved issues as w ell as to capitalize on our ability 
to better control the electron beam ejection process. 
High time resolution observations of plasma waves 
and particles and potential oscillations of beam-emit- 
ting platforms are needed to understand the estab- 
lishment of the payload potential and the sheath struc- 
ture loni/ation processes in the sheath region due to 
out gas. thruster emissions or hollow cathode oper- 
ations also need to be studied. Calculations indicate 
that instabilities may be induced by the ionization 
(Cook and Katz. 1988). 

One of the natural processes that will need more 
studv in the future regards the generation of elec- 
tromagnetic radiation from electron beams. While the 
consensus is that radiation observed in past exper- 
iments is a result of Cherenkov radiation, it is still 
debated how the radiation current structure is formed. 
It is important to answer this question as it pertains 
to the generation of VLE waves in the ionosphere 
from natural processes. The VLE whistler mode is the 
most dominant wave mode observed in the mag- 
netosphere. and these waves control the energetic elec- 
tron population (Kennel and Petschek. 1966). 

The ideas and experiences of the past will be put to 
practice in the .47/..IS-1 space shuttle mission (1992) 
which will fly the SEP AC electron beam experiment 
<7 kV. 1,6mA). For the first time in ionospheric beam 
experiments, a hollow cathode will be used to control 
the spacecraft potential. Thus the beam should escape 
with almost its complete and known energy. Optical 
emissions from the natural aurora as well as from the 
SEP AC electron beam, as it ionizes the Earth's upper 
atmosphere, will be observed by the sensitive on- 
board cameras and spectrographs of the AEPI exper- 
iment. The dynamics of the ionization structures cre- 
ated in the upper atmosphere by the beam will be 
studied by ground-based radars, and attempts will be 
made to observe VLF radiation on the ground. (The 
beam can be modulated up to 5 kHz.) 
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Jhc structure of the radiating beam electron current 
can possibly be deduced from detailed observations 
of the electromagnetic Held in the vicinity of the beam. 
Such measurements are planned lor the ( HARith- 2b 
sounding rocket experiment (1992). The elec- 
tromagnetic radiation from the beam will be observed 
from a sub-payload equipped with a sophisticated 
wave experiment measuring simultaneously three elec- 
tric and three magnetic field components. The exper- 
iment will also allow the study of the limits to vehicle 
neutralization by thruster emissions. 

Finally we point out that the recent development ot 
comparatively small electron linear accelerators now 
makes possible a whole new class of ionospheric modi- 
fication experiments (Banks el ai. 1991). With elec- 
tron beams in the megaelcctronvolt range it is possible 
to create substantial ionization from the source alti- 
tude down to an altitude of 40 km. in the region from 
40 to 60 km. the modification of the fair weather 
electric potential may result in electric fields of a 
sufficient magnitude to cause breakdown and dis- 
charge in the beam column. The result is the initiation 
of lightning discharges along the initial beam column 
and associated VLF wave generation. 
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The dynamics of low-/? plasma clouds moving perpendicular to an ambient magnetic field m 
vacuum and in a background plasma is simulated by means of a three-dimensional, electromag- 
netic, and relativistic particle simulation code. The simulations show the formation of the 
charge sheaths at the sides of the cloud with the associated polarization electric field which facili- 
tate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused 
by the larger ion Larmor radius, which allows ions to move ahead and Lag behind the electrons as 
they gyrate. Results on the cloud dynamics and electromagnetic radiation include the following: 
(1) In a background plasma, electron and ion sheaths expand along the magnetic field at the same 
rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. (2) Sheath 
electrons are accelerated up to relativistic energies. This result indicates that artificial plasma 
clouds released in the ionosphere or magnetosphere may generate optical emissions (aurora) as 
energetic sheath electrons scatter in the upper atmosphere. (3) The expansion of the electron 
sheaths is analogous to the ejection of high- intensity electron beams from spacecraft: a stagnation 
point is reached where the expansion is halted, allowing only a fraction of the electrons to escape. 
These electrons may have energies exceeding the initial beam energy (sheath electrostatic poten- 
tial energy). (4) Second-oider and higher-order sheaths are formed which extend out into the 
ambient plasma These sheaths are curved, the curvature increasing with the order of the sheath. 
(5) The formation of the sheaths and the polarization field reduces the forward momentum of the 
cloud. Furthermore, as the cloud moves across the field, sheath particles are peeled off and left 
behind, while fresh particles continue to drift into the sheaths in order to maintain these. This 
process requires oontinued energy supply which further decelerates the cloud. (6) The coherent 
component of the particle gyromotion (the particles are bom with a coherent component at t— ■ 0) 
is damped in time as the particles establish a forward directed drift velocity. In addition, P*rt*ae» 
undergo orbital phase mixing associated with the establishment of the polarization electnc field. 
For a dense cloud in vacuum the damping occurs very quickly, while for clouds in a background 
plasma the polarization fields are shorted, leaving the particles to gyrate as single particles. In this 
case the coherent component is almost undamped. (7) The coherent particle gyrations generate 
electromagnetic radiation. The simulations indicate that artificial plasma clouds released m space 
will radiate clow to the cloud electron or ion cyclotron /upper hybrid frequencies. 


1. Introduction 

The interaction of plasma clouds, or streams of plasma, 
with ambient magnetized plasmas is a basic process found 
in areas ranging from astrophysics to laboratory plasma 
physics. Examples of such interactions include galactic 
plasma jets [ Yusef-Zudeh et aL, 1984; Arp, 1986; Perley et 
at, 1984], the expansion of cometary plasmas into the solar 
wind [Haerendel, 1987], the dynamics of artificial comets 
[Mende et at, 1989; Baker et aL, 1989; Winske, 1989], 
and heating of tokamak plasmas by the injection of plasma 
beams [Wickman and Robertson, 1983]. Plasma cloud dy- 
namics are also important in experiments studying the criti- 
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cal ionization velocity (CTV) phenomenon [Alfv^n, 1954; Pa- 
padopouloi, 1985; Torbert, 1990J recently performed in the 
Earth's ionosphere and magnetosphere as part of the CRIT 
1 program [Brenninj et aL, 1991] and the CRRES program 
[Reasoner, 1992]. Additional areas of cloud physics include 
the phenomenon of outgassing from large, manned space 
structures. Observations indicate that a large cont am i n a n t, 
partly ionized, water vapor cloud travels with the space shut- 
tle [Kurth and Prank, 1990]. The ionization of the water 
vapor comes about througi charge exchange with ambient 
0 + . The cloud dimension is many HjO + Larmor radii per- 
pendicular to the shuttle velocity vector and extends at least 
300 m behind the shuttle [Patenon and Prank, 1989]. For 
many experiments the outgased plasma cloud represents a 
serious contamination. Thus one of the issues raised by ob- 
servations is to what extent the plasma cloud moves with 
the shuttle [Katz et aL, 1984; Eccles et aL, 1989]. 
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The examples listed above cover a wide range of plasma 
parameters, and consequently theories have been developed 
in the past for various parameter regimes. It has been found 
that plasma beams can propagate in the case where the 
transverse beam width is wide compared to the beam ion 
Larmor radius and fi >> 1 {fi is the ratio of the beam pres- 
sure to the magnetic field pressure) [ Chapman , I960]. This 
regime can be described by MHD theory. The case of high-^ 
narrow beams also allows some cross- field propagation ac- 
cording to hybrid simulations [Papadopoulos et a l. t 1988]. 
However, space and laboratory experiments concur in find- 
ing that the magnetic diffusion of the magnetic field back 
into the beam is much above its classical value. It has been 
suggested that the high diffusion rate is caused by anoma- 
lous resistivity associated with electrostatic wave instabili- 
ties [Mishin et al., 1986]. Finally low-fi beams have been 
shown to cross- field propagate for narrow beams provided 
their density is above some threshold [Peter and Rostocker , 
1982]. 

In this paper we study the propagation of a low-fi cloud 
with a transverse dimension that is larger than the ion Lar- 
mor radius. When a low -fi cloud of plasma is made to stream 
perpendicular to a magnetic field Bo with a bulk velocity tfo, 
one's first guess is that the field will arrest the cloud. How- 
ever, ions and electrons will attempt to rotate in opposite 
directions, with the result that charge layers are formed at 
the sides of the cloud. The electric polarization field Bp 
created inside the cloud by these lateral charge layers is ori- 
ented such that the Bp x Bo drift of the cloud particles is 
in the same direction as the initial bulk velocity. 

A simplified condition for cross-field propagation can be 
found by noting that the crces-field drift velocity v\ induced 
by the electric polarization field is given by 

v\ ~ Bp x B 0 / Bo (I) 

The magnitude of the electric field depends on the charge 
accumulated in the sheaths on the sides of the cloud. For an 
infinite slab of plasma, an estimate of the polarization field 
can be found from Poisson’s equation 

Ep = en c ALy/co (2) 

where c is the elementary charge, n c is the cloud particle 
density, and A Ly is the thickness of the charge layers. As- 
suming the maximum possible sheath thickness to be one ion 
Larmor radius r^,-, a condition for cross-field propagation is 
found from (1) and (2) as 


where u/p, is the cloud ion plasma frequency and u> c ci the 
cloud ion cyclotron frequency. The parameter g* is also equal 
to c/vx where c is the velocity of light and v a the Alfv^n 
velocity. As will be shown later, a one- dimensional cloud 
(an infinite slab of plasma) with g, = 1 is able to cross-field 
propagate with half the initial velocity. Large values of g, 
allow even a three-dimensional cloud to maintain its motion 
across B 0 with almost the full initial velocity. 

The stability of the charge layers has been studied ana- 
lytically and by means of one- and two-dimensional parti- 


cle simulations by Cai et al. [1990] and references therein. 
The cross-field propagation and electromagnetic radiation 
from plasma beams have been studied analytically [Schmidt, 
1960; Borovsky , 1987; Roussel-Duprt and Miller , 1992a, 6], 
while two-dimensional, electrostatic particle simulations 
have been undertaken in order to determine particle diffu- 
sion [Okudo and Hiroe , 1987] and beam dynamics [Livesey 
and Pritchett, 1989]. The dynamics of plasma clouds has 
been explored by one-dimensional electrostatic particle sim- 
ulations [ Galvez , 1987] and two-dimensional electrostatic 
simulations [Galvez et al., 1989; Galvez and Borovsky, 1992]. 
Results from these studies include the rapid expansion of the 
polarization charge layers formed at the sides of the cloud 
along the magnetic field, which tend to deplete the charge 
layers and inhibit the motion of the cloud across the mag- 
netic field, as well as the erosion of the front of the cloud and 
the peeling off of the charge layers as the cloud propagates 
across the field. 

It is clear, however, that the problem inherently is three- 
dimensional and that the presence of a background plasma 
dramatically modifies the dynamics. While much of the 
physics is electrostatic in nature, internal currents as well as 
external currents between the cloud charge layers and the 
ambient plasma generate electromagnetic radiation. The 
radiation from the external currents is analogous to the ra- 
diation from large conductors in space where the cross-field 
motion induces a potential difference between the two ends 
of the structure which, in turn, collect currents from the 
ambient plasma [Donohue et al, 1991]. Electromagnetic 
radiation from plasma clouds is important from an observa- 
tional point of view in the sense that remote observations 
of wave emissions can serve as a diagnostic tool probing the 
cloud dynamics provided the radiation properties of clouds 
are understood. This argument applies both to artificial and 
natural occurring plasma clouds. 

We have therefore extended the simulations referenced 
above by performing fully three-dimensional, electromag- 
netic, and relativistic particle simulations of clouds prop- 
agating both in vacuum and in a background plasma. This 
allows the simultaneous study of the loss of cloud momentum 
in the direction of propagation (x direction), the formation 
of charge layers at the sides of the cloud with the associated 
polarization electric field (y direction), and the expansion 
of the charge layers along the magnetic field ( z direction). 
In addition to the particle dynamics, we also explore the 
radiation properties of clouds. For this first study, a sim- 
ple cloud geometry and initial conditions have been chosen, 
namely a rectangular box of uniform density which is given 
a push across the magnetic field at t = 0. This choice has 
been made in order to better understand some of the funda- 
mentals of cloud dynamics and to evaluate our results rel- 
ative to one- and two-dimensional electrostatic simulations 
using a similar geometry and initial conditions [Galvez, 1987; 
Galvez et ol, 1989; Galvez and Borovsky , 1992]. In section 
2 we present the simulation results while section 3 contains 
a discussion. 

2. Simulation Results 
2.1 Simulation Parameters 

The simulation code used in this study is a further devel- 
opment of the TRISTAN code [Buneman et al., 1980; Pevatt 
et al 1980] originally created for the CRAY-1 and CRAY- 
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XMP. Several new features have been implemented in order 
to increase the speed and versatility of the code. These 
include the elimination of Poisson’s equation and its time- 
consuming solution. Instead, the fields are directly updated 
through current deposition of the particles [ Villasenor and 
Buneman, 1992]. Furthermore, radiating boundary condi- 
tions for the fields are introduced at all boundaries using 
a first-order Lin dm an approximation [Lindman, 1975], and 
finally, Fourier transforms are completely eliminated by pro- 
cessing and filtering the fields locally. The particles have pe- 
riodic boundary conditions. These are not consistent with 
the field boundary conditions; however, as will be clear later, 
this is a problem mainly when the cloud sheath electrons hit 
the boundaries or when the ambient plasma currents to the 
sheaths are exhausted. When the boundary conditions in- 
terfere with the physics, the simulations are stopped. 

The simulations are performed on the CRAY-2 located at 
the National Center for Supercomputing Applications, Uni- 
versity of Illinois. The simulation domain with the plasma 
cloud is shown in Figure 1. The domain is a mesh of si ze 
92 A x 46 A x 185 A where A is the grid unit. The cloud 
dimension is 40 A x 12 A x 12 A. The magnetic field is in the 
z direction, and at f=0 the cloud is pushed across the field 
with a velocity v 0 in the x direction. 

Four cases are simulated, corresponding to two sets of 
cloud and background plasma parameters. The background 
is either a vacuum or a plasma. For the cloud we have chosen 
two plasma densities corresponding to — 1.25 and 2.50. 
The densities relate to the background plasma density as 
n c fn b = 2 and 8, respectively. The cloud electron thermal 
velocity normalized to the velocity of light, equals 

0.02 for the rarefied cloud and 0.04 for the dense cloud. 
This choice ensures that the cloud Debye length is 0.25 A 
for both clouds. The cloud ions have the same mass as the 
background plasma ions, and the ion to electron mass ratio is 
rm/m e = 16. Ions and electrons have equal thermal energy. 
For the background plasma, v\ h (c = 0.02, where is the 
electron thermal velocity. The cloud is given an initial bulk 
velocity v 0 at t=0 in addition to the thermal velocity. The 
magnitude of this velocity is v c /v b h = 1.6. From the above 



Fig. 1. The simulation domain. 


parameters it can be seen that the rarefied, cold cloud has 
a drift velocity in excess of the thermal velocity while the 
reverse is true for the dense, hot cloud. The ratio of the 
particle pressure (thermal ^dynamic) to the magnetic field 
pressure 0 is 2.2 x 10 -3 for the rarefied cloud and 1.6 x 10 
for the dense cloud. The dimension of the cloud is such that 
the ion L armor radius due to the drift velocity is = 
2A. The above parameters lead to the following values of 
characteristic plasma periods: the ion cyclotron period T& 
= 785 time steps, the electron cyclotron period Tee = 49, 
and the background electron plasma period T$e = 222. For 
the rarefied cloud we have = 628 and 1f h = 490, where 
is the ion plasma period and is the ion upper hybrid 
period, while for the dense cloud, T £ =314 and TJ^ — 233. 
The ion upper hybrid per iod is defi ned from the ion upper 
hybrid frequency u>f h = 

2. 2 Sheath Formation 

The charge sheaths as they have formed at t = Td/ 4 are 
shown in Plate 1 for the case of qi = 1.25 and vacuum. The 
top panel shows the space charge in the center four layers in 
the xy plane. Negative charge is blue, and positive charge is 
red. The arrows indicate the direction of the electric field in 
that plane. The lengths of the arrows have been scaled in a 
way which enhances the weaker fields in order to see the field 
structure at some distance from the cloud. The background 
magnetic field is pointing out of the plane, and the cloud is 
moving toward higher x values. As the ions have gyrated 
clockwise and electrons counterclockwise, a positive charge 
layer is created on the right side of the cloud (right side as 
defined relative to the direction of motion of the cloud, i.e., 
low y values), and a negative layer on the left side (high y 
values). In addition, & positive layer has formed at the head 
of the cloud, and a negative layer at the tail. These layers 
are a result of the larger ion L armor radius, which allows 
ions to move ahead of electrons across the magnetic field. 
The electric field inside the cloud is pointing predominantly 
in the y direction, while the field outside of the cloud has 
a twisted, dipolarlike structure arising from the L-shaped 
charge layers. Inside the cloud some smaller modulation of 
the space charge density has formed. 

The charge density in the xz plane is shown in the bot- 
tom left panel. The cloud is moving horizontally and in the 
positive x direction, and the magnetic field is vertical and in 
the z direction. At the front of the cloud is seen the positive 
ion sheath developed by ions moving ahead of electrons, and 
in the tail the negative electron sheath resulting from elec- 
trons lagging behind ions. The electron sheath has at this 
point started to expand along the ambient magnetic field as 
a result of the acceleration by the parallel component of the 
electric field. Because of the higher ion mass, the ion sheath 
has expanded much less. 

Close to the cloud the electric field structure is dipolar in 
nature, while at larger distances the field changes character 
and periodically changes orientation. As we will show later, 
the oscillating fields are those of an outward propagating, 
right-hand electromagnetic wave generated by the current 
induced from the coherent motion of electrons in the cloud. 
A dipole electric field decreases with distance r as r~ 3 while 
radiation fields decrease as r~ 2 . Therefore the radiation 
fields dominate over the electrostatic field at some distance 
from the cloud. 

The charge density in the yz plane is shown in the bottom 
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Fig. 2. The location of cloud iona projected onto the xy f xz t and yz planes at t = T« foe the dense cloud («= 2.50) 
in vacuum. 


right panel. In this perspective, the dond is moving out of sheath is beginning to peel off at the tail of the dead. The 
the plane with the magnetic field vertical and in the * direc- polarisation electric field decreases through the sheaths from 
tion. The electric field is approximately dipolar outside of its full value at the inner edge to zero at the outer edge of the 
the doud. As was the case for the electron tail sheath, elec- sheaths. Therefore the average drift velocity of the particles 
trons have begun their expansion along the magnetic field, in the sheaths decrease from the bulk velodty at the inner 
which elongates the negative charge layer in the z direction, edge to zero at the outer edge [Botvvtky, 1987; Cat tt al.. 
Note that the strongest parallel electric fields are found at 1990]. Notice, also, that a region of ion density depletion has 
the tips of the charge layers. Thus dectrons located here developed which is dose to and aligned with the primary 
will experience the most rapid acceleration. ion sheath. As we shall see later, secondary space charge 

The evolution of the dense cloud in vacuum at a time t = sheaths develop within the cloud and, in the case of a finite 
T ei is illustrated in Figure 2, which shows the ion positions background plasma density, also in the background plasma, 
projected onto the xy plane (top), the xz plane (bottom The expansion of the ion sheaths along the ambient mag- 
left) and the yz plane (bottom right). As seen from the netic field is shown in the bottom frames. In the bottom 
top frame, the cloud has broadened at the head, and the ion left frame the expansion at the tail end of the doud is par- 
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ticularly noticeable. As will be shown later, the bulk of the 
ions oscillate approximately with the ion upper hybrid pe- 
riod. From the values of the characteristic time constants 
given earlier it is seen that the ions at this time have ex- 
perienced more than three oscillations (for the dense cloud, 
Ti h IT c i = 0.30). Ions have therefore been ahead and be- 
hind electrons several times, and thus positive and negative 
charge layers have been formed alternately at the head and 
tail of the cloud. The expanding ion layer seen in the tail 
is at a location of the rear edge of the bulk of the cloud 
and corresponds to the narrow layer of ion density enhance- 
ment seen in the top frame. The ions that have peeled off 
the ion sheath are seen as a diffuse distribution behind the 
expanding tail layer. 

The expansion of ions in the ion sheath at the side of 


the cloud is seen in the bottom right panel. The ion layer 
has a small dispersion, with ions far from the cloud located 
at y values closer to the center y coordinate of the cloud. 
These ions were the first to leave the cloud in their expansion 
along the magnetic field, and their location in y indicates the 
position where the charge layer first forms. With time, the 
cloud has broadened, in particular in the front portion. 

The electron spatial distribution is shown in Figure 3. In 
the top panel is seen the peeling off of electrons in much 
the same way as was observed for the ions. Most notice- 
able, however, is the expansion of the electrons in the nega- 
tive electron charge layer out to the edge of the simulation 
domain shown in the bottom two panels. In the process 
the most energetic electrons reach relativistic energies. The 
expansion of the lateral electron charge layer is in qualita- 


ELECTRONS 




Fig. 3. The location of cloud electrons projected onto the xy, and yz planes at t = T* for the dense cloud in 
vacuum. 
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tive agreement with two-dimensional, electrostatic simula- 
tions [Galvez and Borovsky , 1992]. Notice also (bottom right 
frame) that electrons have accumulated in the ion layer in 
their attempt to neutralize the ion positive space charge. 
Furthermore, a secondary electron layer has formed close to 
the center of the cloud. The interior of the cloud is then 
expected to be quite dynamic with large-amplitude density 
and space charge oscillations. 

The bottom left frame of Figure 3 gives a side view of the 
cloud and the electron sheath. The sheath exhibits a small 
dispersion, the outermost electrons lagging behind the cloud. 
The result is the formation of an electron “wing” slanted 
with an amount corresponding to the cross-field motion of 
the cloud and the field-aligned motion of the electrons. 

Results from two-dimensional simulations (our xy plane) 
of clouds with qf = 0.5, 10, and 32 at a time around one ion 
cyclotron period are reported by Galvez et al. [1989]. These 
results are comparable to those shown in Figures 2 and 3 in 
the sense that they exhibit a broadening of the pulse front 
and the peeling off of the charge layers. Furthermore, Galvez 
et al. find the development of vortices at the front end of 
the cloud at the edge of the electron sheath, and vortices in 
the tail end of the cloud. The front vortex was found to be 
charge neutral and a result of a bulk plasma drift at the head 
of the cloud in a direction toward the main negative charge 
layer. This motion is an £ x S drift due to the electric field 
component antiparallel to the direction of motion generated 
by the positive space charge at the head of the cloud. The 
vortices at the tail of the cloud were not charge neutral but 
a result of the peeling process. In addition, for the high- 
density cloud they found the development of the diocotron 
instability [Buneman, 1961] in the electron sheath. 

The tendency to develop these vortices and instabilities 
as studied by Galvez et al. [1989] increased with plasma 
density. In our own simulations, with qf — 6.25, these phe- 
nomena were suggested, but little developed at t = T ei . 
We propose that this is because our simulations are three- 
dimensional, which allow for a more rapid decrease of the 
electrostatic fields with distance. Furthermore, since the 
simulations by Galvez et al. [1989] were performed in a 
plane perpendicular to the magnetic field, electrons were 
not allowed to expand out of the charge layers. 

For cloud propagation in a background plasma the plasma 
currents are able to short the space charge sheaths for both 
the rarefied and the dense cloud. Thus space charge layers 
are formed to a lesser extent, and the cross-field motion is 
inhibited. As a result, the features seen in Figures 2 and 3 
are much less developed. 

The charge distribution in the yz plane at t=T c i is shown 
in Plate 2, left panel, for the rarefied cloud in vacuum, The 
charge layers at the sides of the cloud have expanded to 
the extent indicated by the spatial particle distributions. 
The ion sheath has expanded very little, while the electron 
sheath has almost reached the boundary of the simulation 
domain. The electric field structure has now departed from 
the dipolarlike structure seen in Plate 1 at t = T c {j 4. The 
field orientation is such that electrons in the electron sheath 
located at y values between about 30 A and 150 A now are 
decelerated, while electrons further from the cloud continue 
to be accelerated. While the longer time duration dynamics 
of the sheaths cannot be studied with this simulation be- 
cause of limitations on the size of the simulation domain, 
we suggest nevertheless that the outermost electrons may 
escape the cloud while some of the electrons closer to the 


cloud eventually will return to the sheath. The situation is 
similar to the case of electron beam emission from space- 
craft. Simulations find here that for beam currents that are 
large compared to the plasma currents, a beam stagnation 
point will develop, allowing some electrons to escape with 
energies in excess of the beam accelerator energy, while the 
bulk of the electrons are returned to the spacecraft [Winglee 
and Pritchett, 1988; Pritchett, 1991]. 

When a background plasma is present, currents are drawn 
to the sheaths, thereby reducing the polarization field. In 
Plate 2, right panel, is shown the charge distribution for 
the rarefied cloud in a background plasma. The two main 
sheaths are present as for the cloud in vacuum. However, the 
electron sheath is less developed, and the ion sheath has now 
expanded to the same degree as the electron sheath. The 
increased ion sheath is a result of electrons being collected 
from the background plasma to the sheath. This electron 
current leaves behind a positive space charge which has the 
effect of elongating the ion sheath. In addition to the two 
primary sheaths, several secondary electron and ion sheaths 
develop. The outermost of these seen on Plate 2 are nega- 
tive electron sheaths which are curved rather than aligned 
with the magnetic field. In fact, the higher the order of the 
sheath, the more curved it is. For instance, next to the pri- 
mary electron sheath, and at higher y values, is a slightly 
curved ion sheath. We may for convenience call this the 
second-order sheath. Next to this is found the third-order 
electron sheath which is highly curved, and similarly for the 
primary ion sheath, which has a slightly curved (in the op- 
posite sense) second-order electron sheath. The third-order 
ion sheath is not visible. 

The higher-order sheaths are developed by electrostatic 
field component is not at a y value corresponding to the loca- 
tion of the charge layers, but at y values which are somewhat 
further from the cloud. As a consequence, ambient plasma 
located further from the sheath is accelerated the most, and 
the curvature mentioned above is created, with a higher de- 
gree of curvature for higher-order sheaths. As an example, 
consider the field from an electric dipole oriented with its 
axis l perpendicular to an ambient magnetic field S 0 Let r 
be the vector from the center of the dipole to the point of 
observation and 9 the angle (f, ?). It can be shown that at 
a given distance r from the dipole center, the electric field 
component parallel to S 0 maximizes for 6 = 45°. If the 
dipole is immersed in an ambient plasma, the charged parti-^ 
cles located along the line B = 45° in the plane defined by l 
and So will experience the strongest acceleration along B 0 . 
As these particles respond to the fields, a new curved layer 
is formed. The fields from this layer in turn create an even 
more curved layer, and so on. 

The cloud electrons in the negative charge layers are ac- 
celerated to relativistic energies. The strong acceleration is 
caused by the sudden push of the cloud across Bo at t — 0 
and the associated impulsive formation of the charge lay- 
ers and electric fields. An estimate of the maximum energy 
of electrons accelerated along the magnetic field due to the 
electrostatic forces in the electron sheath can be found by 
noting that the maximum potential energy of an electron in 
the sheath is 


Wpot = - tV 0 B 0 Ly 


( 4 ) 


where Ly is the dimension of the cloud in the direction along 
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Plate 2. The charge density and electric field in the center four layers of the y z plane at t = T ci for the rarefied 
cloud (<jj=1.25) in vacuum (left) and in a background plasma (right). Positive space charge is red, and negative 

blue. 
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the polarization field. Equating the kinetic energy of an plied from the cloud kinetic energy. As a consequence the 
electron accelerated out of the sheath with the potential cloud slows down from its initial velocity v a . The loss rate 
energy in the sheath, the maximum kinetic energy becomes of particle energy density can be found from 


U! max ^ TTljLy 
w% m e r Ll 

where Wq = The acceleration out of the sheath 

results according to (5) in electron energies which are equiv- 
alent to the ion bulk kinetic energy multiplied by a geomet- 
ric factor. For clouds which are wide compared to the ion 
Larmor radius, very large energies can be achieved. 

The maximum eneigy given by (5) is 89.1 for the dense 
cloud and 75.0 for the rarefied cloud. These values are of the 
same order of magnitude as found in the simulations. The 
energy of the fastest electron relative to the initial drift en- 
ergy was for the dense cloud 80.7 and 241.3 in a background 
plasma and vacuum, respectively, while it was 21 .9 and 100.9 
for the rarefied cloud. They are lower than the vacuum cases 
and higher than the background plasma cases. When deriv- 
ing (5), the expansion of the electron charge layer and the 
continuous supply of charge from the cloud, which supply a 
constant push, were not accounted for. Again we point out 
the analogy with electron beam experiments. Electrons at 
forces from lower-order sheaths. As pointed out by Galvez 
and Borovsky [1992], the maximum of the parallel electric 
the beam front are accelerated by electrostatic forces to en- 
ergies which exceed the initial accelerator energy. Therefore 
the maximum energy in the vacuum case is expected to be in 
excess of the predictions of (5). With a background plasma, 
the fields are partly shorted by ambient plasma currents, 
and the acceleration is less. 

2.3 Cloud Energy and Momentum 

The generation of electric and magnetic fields and the 
energization of sheath particles require energy which is sup- 


d -K = ].£-R (6) 

dt 

where J is the current density and £ the electric field. The 
energy loss due to electromagnetic radiation from the cloud, 
R, is relatively small compared to the energy loss incurred in 
establishing the charge sheaths and will be discussed later. 

The current required to supply the charge layers at the 
sides of the cloud and the associated polarization electric 
field are dominant in this energy transfer. The polarization 
electric field is oriented in the positive y direction while the 
current supplying the charge layers at the sides is oriented 
in the negative y direction. The process of forming and 
sustaining these layers therefore reduces the cloud kinetic 
energy. The charge layers at the head and tail of the cloud 
alternate in polarity during an ion oscillation. The current 
in the x direction supplying these layers is 90° out of phase 
with the field, and the energy will therefore transfer both 
from the cloud kinetic energy to the fields as well as back 
again. This, together with the smaller area of the head and 
tail surfaces of the cloud relative to the sides, results in a 
comparatively small effect on the cloud energy budget from 
these layers. 

The coherent components of the cloud ion and electron 
velocities are found by averaging the velocities of all cloud 
ions and electrons. The x components of the coherent veloc- 
ities are shown as a function of time in Figure 4, left frame, 
for the rarefied cloud in vacuum. In the figure, the velocities 
have been normalized with the initial drift velocity v 0 in the 
x direction. The ion velocity (dashed curve) oscillates with a 
period between the ion upper hybrid period (Tf^/Td — 0.62) 
and the ion cyclotron period. The oscillations are damped in 



Fig. 4. The average cloud ion velocity (dashed line) and cloud electron velocity (solid line) in the x direction as 
a function of time for the rarefied cloud in vacuum (left) and in a background plasma (right). The velocities are 
normalized to the initial drift velocity. 
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time, and the average velocity approaches 0.2 at the end of 
the run. The electrons oscillate at the electron cyclotron fre- 
quency, these oscillations also being damped in time. Note 
that since w C e > we cannot distinguish between the 
electron cyclotron period and the electron upper hybrid pe- 
riod. The oscillations in the velocities are damped partly 
because individual ions and electrons gyrate less as the po- 
larization field is getting established, and partly because the 
particles lose phase coherence. 

The characteristics seen in Figure 4, left frame, are to be 
expected from a cloud that is marginally able to cross-field 
propagate as defined by (3) and are in qualitative agreement 
with two-dimensional simulations performed in the past. A 
cloud as defined initially in our simulations will pass from 
an initial stage where it is given a push across B 0 > through 
a transitory stage where the charge layers are being estab- 
lished and the cloud particles undergo significant coherent 
gyrations, to a relaxed stage where the layers are established 
and the cloud ions and electrons drift with a common (lower) 
velocity across B 0 with very little or no coherent gyrations. 
For a cloud that is marginally able to cross-field propagate, 
the bulk velocity in the final stage is much smaller than the 
initial velocity. As mentioned in the introduction, a one- 
dimensional cloud (slab) will in the latter stage cross-field 
propagate with a velocity that is half the initial velocity if 
qi = 1. For the case shown in Figure 4, left frame, g, = 1.25, 
which in the one-dimensional case should result in a final 
drift velocity slightly above 0.5 v 0 . However, in three di- 
mensions, the charge layers are allowed to expand, which, 
as we will show, leads not only to a smaller final drift veloc- 
ity, but also to a steady, slow decrease of this velocity. 

The coherent velocities for the rarefied cloud in a back- 
ground plasma are shown in Figure 4, right frame. After an 
initial damping, electrons continue undamped oscillations 
at the electron cyclotron frequency with very little net drift 
across the field. Ions oscillate at first with the ion cyclotron 
period, and at about t = 1.4 T c i, ions begin to oscillate with 


a period closer to T* h . The reduction of the polarization 
electric field inside the cloud due to the currents from the 
ambient plasma to the charge layers makes the cloud par- 
ticles oscillate at the cyclotron frequencies much as single 
particles. In other words, the shorting of the polarization 
electric field reduces collective effects. The oscillation pe- 
riod of the ions changes from the cyclotron period to the 
upper hybrid period because the ambient plasma currents 
have been exhausted. These currents flow in thin sheaths 
along the magnetic field, and at the time of change in pe- 
riodicity, little charge is left because of the finite extent of 
the simulation domain. At this point the cloud ions start to 
oscillate as if in vacuum. The effect of the boundaries sets 
the limit in time that the simulations can be run. 

The velocity for the dense cloud in vacuum is shown in 
Figure 5, left frame. The electron oscillations at the electron 
cyclotron frequency are strongly damped and disappear af- 
ter a few oscillations. The low-frequency oscillation in the 
ion and electron velocities is now close to the ion hybrid 
period ( Ti h /T ei = 0.3), and the loss of forward momentum 
is less than for the rarefied cloud. After an initial decrease 
within the first half of an oscillation period, the ion veloc- 
ity oscillation is slowly damped while oscillating around a 
slowly decreasing average value which reaches about 0.4 at 
the end of the run. 

Finally, the velocity for the dense cloud in a background 
plasma is shown in Figure 5, right panel. The ions now 
oscillate much slower and closer to the cyclotron period. 
At about t = 1.1 T a - the ambient plasma currents are ex- 
hausted, and the oscillation period decreases, as was the case 
for the rarefied cloud. After an initial damping, electrons 
continue their almost undamped oscillations at the electron 
cyclotron frequency. The oscillations are of a smaller ampli- 
tude than was the case for the rarefied cloud. The simulation 
is terminated shortly after the influence of the boundaries is 
noted in the cloud velocities. 

Some of the characteristics seen in Figures 4 and 5 can 
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Fie 5 The average cloud ion velocity (dashed line) and cloud electron velocity (solid fine) in the * direction as 
function of Sfor the dense cloud in vacuum (left) and in a background plasma (right). The velocities arc 


normalized to the initial drift velocity. 
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be understood from considerations of the energy budget of 
the clouds. To simplify the situation we estimate the en- 
ergy flow associated with the formation and maintenance of 
the primary charge layers at the sides of a cloud in vacuum. 
First we note that the electron charge layer rapidly expands. 
While some of the first and fastest electrons accelerated out 
of the sheath may be able to escape the cloud completely, 
the electric field structure as shown in Plate 2 seems to indi- 
cate that there may be a stagnation point in the expansion 
where £|| = 0. It is not known if the stagnation point is 
stationary or if it oscillates by contracting and expanding 
in time. However, the slow decrease of ion momentum after 
the initial drop mentioned above seems to indicate a steady 
loss of momentum. In the following we then assume that 
the layer reaches some quasi-steady, average position in the 
direction along B 0 . 

Let the original cloud dimension be L X} L y ,L z in the 
x, y, z directions, and let the dimension of the extended elec- 
tron layer along the z direction be M z < As this layer is being 
formed, the charge required to increase the thickness of the 
layer by dL y is dQ = neL x dL y M z , where n is the electron 
density in the layer. This charge can be expressed in terms 
of the associated increment of the polarization electric field 
in the cloud dEp by using (2) and the assumption of an 
infinite slab: 


where q 2 = q\ 1 + q\ = r/v 2 A . 

The expression for vi has been derived under the assump- 
tion that the particles contain little gyrational energy and 
thus constitutes an upper limit for the cross-field drift at a 
time after the establishment of the charge layers. However, 
particles do retain some gyromotion. A second estimate of 
the final cross-field velocity comes from arguments of conser- 
vation of canonical momentum of the cloud particles. From 
Lorentz’ equation one finds that midv x /dt = eB 0 dy/dt. 
The y-ward displacement of the ions is then associated 
with a reduction in velocity in the x direction. The same 
is true for electrons, which are displaced to the opposite 
side. The displacement exposes charge layers of thickness 
dvx(mj 4- m c )/eB 0 on the two sides and results in a polar- 
ization electric field Ep = n c (mj 4 m e )dv x / B 0 c 0 . The final 
drift velocity is then 


or 


dv x n c (mi 4 m e ) 

— zm 


(13) 


Vx 

Vo 


r 41 


(14) 


dQ — tod Ep L x M z 


(?) 


The current density in the cloud required to supply this 
charge flows in the negative y direction and is given by 


T dE' p Mz 

J ± = t °-jr it 


( 8 ) 


The corresponding energy density flow is then 


dW_ _ l dE p 2 
dt ~ V° dt it 


( 9 ) 


where we have used (6). The total energy density required 
to establish the sheaths then becomes 


W = 


( 10 ) 


where Ep is an average, quasi-steady polarization field at the 
time when the cloud has established its main polarization 
layers. 

In the hypothetical case of no expansion of the charge 
layer, M z /L z = 1 and (10) becomes the familiar expression 
for the energy density of an electric field in vacuum. Assum- 
ing that most of the energy flow goes into establishing the 
lateral charge layers, conservation of energy then requires 


Here the influence of the expansion of the electron charge 
layer has been ignored, while the loss of forward momentum 
is more rigorously accounted for. 

We propose now that the initial decrease of momentum 
of the cloud primarily is a result of the energy required to 
push charge into the lateral change layers during the first 
half ion oscillation. We also propose that the later, slower 
deceleration of the cloud is caused by the current flow in the 
cloud required to maintain the layers as the cloud moves 
across the magnetic field. Some energy is also radiated from 
the cloud, in particular in the transition phase, as will be 
discussed below; however, it constitutes a relatively small 
energy loss. From (12) an estimate of the degree of ex- 
pansion of the electron charge layer ( M z / L z ) can be found. 
Although difficult in the case of the rarefied cloud, we vol- 
unteer an estimate of the initial velocity decrease of this 
cloud in vacuum (Figure 4, left frame), namely, v\/v 0 — 0.4, 
which leads to M z /L z ~ 8.1. The extent of the charge layer 
along B 0 is then 97. 2A, which compares favorably with the 
distance between the two stagnation points (£y =0) seen 
in Plate 2. For the dense cloud, v\/v 0 ~ 0.6, which gives 
Mz/Lz ~ 11 . 1 - 

As the cloud moves across the magnetic field, the sheaths 
peel off as seen in Figures 2 and 3. The average current 
density required to maintain the sheaths is given by 



vi M z 
L x L z 


( 15 ) 


i(m; 4 m t )v 2 0 n c ^ i(m; 4 m e )v\n c 4 -e 0 Ep-~- ( 11 ) 

Substituting v\B 0 for E p and rearranging the terms, one 
finds 


The factor £ comes from the fact that the forward velocity 
of the electrons in the sheath decreases roughly linearly from 
a value of v\ at the inner edge of the sheath to 0 at the outer 
edge. The time constant associated with the deceleration of 
the cloud due to the maintenance of the layers becomes 


/ \ 2 2 

~ ? 

\V 0 / q 2 + Mz/Lz 


(12) 


q? L X L Z v 0 ~ 

r = - rz T ci 

X TliMz Dl 


( 16 ) 
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Fig. 6. Relative amplitude of the y component of the electric field along an axis in i the z direction through the 
center of the rarefied cloud in vacuum (left frame) and in a background plasma (right frame). Along the vertical 
axis is amplitude and time, and along the horizontal axis is z. 


where r is defined as r = /(dv\j dt). For the rarefied 

cloud, vi/vo ~ 0.2 and M z /Lz ^ 8, which gives r ~ 6.1 T«. 
For the dense cloud, v\/vo cz 0.5 and Mz/Lz — lli which 
gives t cz 7.3 T c i . These numbers are in reasonable agree- 
ment with the observed deceleration. With our assumption 
of the initial velocity decrease mentioned earlier, both clouds 
decelerate by 0.2 in the remaining time up to 2 T C1 . Assum- 
ing that this second slower decrease begins at a time when 
the ions have completed 180° rotation, this corresponds to 
r = 8.25 T c i the rarefied cloud and r = 8.74T c j for the 

dense cloud. 

2.4 Electromagnetic Radiation 

We now turn to the electromagnetic radiation from the 
clouds. In Figure 6, left frame, is shown the amplitude of 
the y component of the electric field Ey along an axis in 
the z direction through the center of the rarefied cloud in 
vacuum. The vertical axis represents both relative ampli- 
tude and time. In the center of the frame is seen the growth 
and later decay of the polarization electric field inside the 
cloud. Electromagnetic radiation at the electron cyclotron 
frequency is also generated, which propagates into the vac- 
uum from the cloud with the velocity of light. These waves 
are generated fairly strongly in the initial phase but are 
damped out relatively quickly. In Figure 6, right frame, 
is shown E y for the rarefied cloud in a background plasma. 
The amplitude is here again relative, but scaled in the same 
manner in order to be directly comparable to the values 
shown in the left frame. The polarization field inside the 
cloud is now much more structured and with time extending 
out into the ambient plasma. Recall that the initial cloud 
dimension L z is 12A. These structures may be associated 
with the higher-order sheaths which are formed in the cloud 


as well as in the plasma, as mentioned earlier. The elec- 
tromagnetic radiation continues to be generated throughout 
the time shown on the frame in contrast to the case of the 
cloud in vacuum. 

We propose that the radiation from the cloud is gener- 
ated by the ion and electron currents flowing in the cloud, 
which act as virtual antennas. The 1 component of the cur- 
rents can be found from the average cloud particle velocities 
shown in Figures 4 and 5. For the cloud in vacuum, the 
average electron velocity oscillates strongly at the electron 
cyclotron frequency early in the run, but is later damped 
out, while for the cloud in a plasma, the electrons continue 
to oscillate coherently. This is consistent with the high- 
frequency radiation seen in Figure 6, which is damped in 
time for the cloud in vacuum and continues to be gener- 
ated for the cloud in a background plasma. We note at 
this point that the initial conditions imposed on the particle 
velocities are important for the level of radiation. The par- 
ticles are created instantaneously with a common velocity 
across the magnetic field; thus the cloud particles are born 
with a large coherent velocity component. This issue will be 
discussed further in the next section. Electromagnetic ra- 
diation is also generated in the low-frequency regime at the 
ion oscillation period; however, longer simulations in time 
and a different scaling for the display of the field amplitudes 
are needed to see these waves. 

In order to examine the radiation fields further, the ra- 
diation magnetic field was sampled at z = 140A on the 
axis in the z direction passing through the center of the 
cloud. Figure 7 shows a hodogram of the magnetic field in 
the x, y plane for the rarefied cloud in vacuum around a time 
i _ T ci f 4. A value is plotted for every 5 time steps from step 
150 to step 250. The field is rotating in the same sense as 
the electrons and is therefore the field of a right-hand polar- 
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q=1.25; Vacuum 



Fig. 7. Hodogram of the magnetic field amplitude (relative) in 
the xy plane for the rarefied cloud in vacuum. The location is 
center-x, center-y, and 2 = 140 A. A value is plotted every 5 time 
steps from step 150 to step 250. 


ized wave. As noted earlier, the electron oscillation period 
is about 49 time steps. During the 100 time steps shown 
in Figure 7, the field has undergone a little more than two 
complete rotations and thus oscillates at the electron cy- 
clotron frequency. Because of the finite velocity of light the 
source current generating the fields is the current at earlier 
time steps around 50 to 150. As we will show below, the 
coherent electron current in the cloud constitutes the source 
current. Thus the source current of the fields shown in Fig- 
ure 7 is associated with the second and the third electron 
oscillation seen in Figure 4, left frame. 

The radiation from the cloud in vacuum has been esti- 
mated analytically by assuming a source current of the form 

J = J e e iWcet 1 + jf (17) 

where J t is estimated from the coherent motion of the elec- 
trons seen in Figures 4 and 5 (the oscillation amplitude is 
roughly the same in the y direction). The vectors £ and 
jj are unit vectors in the x and y directions. The vector 
potential A is found from 



where the integration is over the cloud volume and k is the 
wave vector associated with the dominant wave oscillation 
frequency w C e \ k is parallel to the position vector r, which is 
the vector from the volume element dV to the point of obser- 
vation. The magnitude of k is given the vacuum dispersion 
relation w C e/k = c. 

The electromagnetic fields are given by 



and the Poynting flux is found from 

3=i.Re[£x#‘] (21) 

As mentioned above and shown in Figure 7, the fields for 
the rarefied cloud in vacuum have been sampled in the simu- 
lation at z — 140A on an axis in the z direction through the 
center of the cloud at t = T c i/4 - At this time, the charge lay- 
ers have not undergone significant expansion, and the dom- 
inant fields are those of the radiation fields generated from 
the currents in the cloud. The vector potential at this time 
has been found by estimating the current from the coherent 
motion of the electrons in the second and third oscillation 
shown in Figure 4, left frame. Assuming the normalized am- 
plitude of the velocity oscillation is 0.38, the vector potential 
has been found by integrating the associated current density 
numerically. The z component of the Poynting vector esti- 
mated from the vector potential has been compared to the 
z component found directly from the sampled fields. The 
agreement is within 10% with S z = 3.27 x 10 7 when esti- 
mated from the vector potential calculation and 3.0 x 10 -7 
from the direct sampling of the fields (in simulation units). 
This good agreement substantiates our hypothesis that the 
fields are generated by the coherent electron current. 

The loss of momentum due to radiation can be estimated 
by integrating the Poynting flux over a surface around the 
cloud. Assuming the average Poynting flux is of the order 
of 1 x 10~ 7 through an area 2 x 92A x 46A, the energy loss 
rate due to radiation is of the order of 1 x 10” 5 (in simula- 
tion units). The time constant of the braking of the cloud 
associated with radiation losses defined as r ro( f = v/(dv/dt) 
is then about 3x 10 3 T c i, which is much longer than the time 
constant found for the maintenance of the charge sheaths. 

Finally we show the dispersion of the radiation from the 
rarefied cloud in Plate 3. The amplitude of the magnetic 
field component B x in decibels is shown as a function of fre- 
quency and wave number k z for vacuum (left frame) and 
a background plasma (right frame). The plot has been 
created by performing a two-dimensional Fourier transform 
(one space and one time) of B x along the 2 axis through the 
center of the cloud for the first 500 time steps. The radiation 
maximizes as expected around with additional lower- 
amplitude radiation generated along a line with a slope cor- 
responding to the velocity of light. The radiation from the 
cloud in a background plasma shown in the right panel is 
also primarily generated close to uf C e , with lower-intensity 
additional radiation along a line with a slope corresponding 
to the velocity of light. In this case, however, radiation is 
damped right at u C t with the result that two branches ap- 
pear. The one coming from low frequencies goes to large k z 
as the frequency approaches u> c e (a resonance), while the sec- 
ond branch starts at low k z values right above u? C e (a cutoff) 
and approaches the slope for the velocity of light for higher 
frequencies. These are the characteristics of the right-hand 
mode in a plasma with u> C e !> The cutoff and resonance 






Plate 3. Wave intensity aa function of w and k, for the rarefied cloud in vacuum (left) and in a background plasma 
(right). 
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are in this case both at o/ Cf . Note that the amplitudes of 
the cyclotron frequency radiation are comparable in the two 
cases although they appear to be of different amplitude due 
to the color coding. This is caused by the difference in dy- 
namic range used in the two plots. 

3. Discussion 

The shapes of the clouds investigated have been chosen 
as simple as possible in order to be able to differentiate be- 
tween the processes and forces acting in different spatial 
directions. Natural or artificial plasma clouds will in gen- 
eral have smoothly varying density through the cloud, and 
the effect of the charge layers at the sides of such clouds will 
be less dramatic. The initial conditions were also chosen 
simple, namely the sudden push at t=0. This condition has 
as a consequence a sudden acceleration of electrons out of 
the sheath and a burst of radiation generated by the phase- 
coherent particles. In this regard the clouds are not in equi- 
librium initially but relax toward an equilibrium state at the 
end of the runs. 

Given our initial conditions the simulations suggest that 
the deceleration of the clouds in vacuum occurs in two steps: 
The first step is a rapid deceleration within the first half of a 
cloud ion hybrid oscillation period caused by the transfer of 
cloud kinetic energy into the formation of the charge layers. 
In the case of the electron charge layer this layer extends 
out to some stagnation distance from the cloud. The second 
stage is the continuous deceleration caused by the current in 
the cloud needed to maintain the charge layers as these are 
peeled off when the cloud moves across the magnetic field. A 
third stage, which includes the expansion of the ion charge 
layers, has not been studied due to limitations in the length 
of a simulation run. 

Associated with the formation of the electron charge lay- 
ers, electrons are initially accelerated to energies which ex- 
ceed the ion drift kinetic energy. In the case of clouds that 
are large compared to the ion Larmor radius, analytical 
considerations indicate that the energy may be consider- 
ably larger than the ion kinetic energy. This acceleration is 
caused by the sudden push of particles across the field. In 
the case of CIV experiments, or other cloud experiments in 
space plasmas, the ionization rate, or cloud formation rate, 
is slow compared to an electron cyclotron period, and the ac- 
celeration is likely to be smaller. Analytical estimates indi- 
cate, though, that if a cloud is wide compared to the ion Lar- 
mor radius and dense enough to cross-field propagate in the 
presence of a background plasma, electrons may be acceler- 
ated to energies which are comparable to or exceed the elec- 
trostatic sheath potential energy. In this respect the elec- 
tron sheath expansion has much in common with electron 
beam ejections from spacecraft. For dense beams a virtual 
cathode will form in the beam at a point where the beam 
electrons stagnate as their velocity away from the spacecraft 
is reduced, due to the electrostatic forces in the beam and 
the spacecraft, which charges to positive potentials [Wxnglee 
and Pritchett , 1988; Pritchett , 1991]. The virtual cathode 
allows only a fraction of electrons to escape the spacecraft; 
however, some of the escaping electrons may have energies 
that exceed the beam accelerator energy. Clouds released in 
the ionosphere or magnetosphere may then directly excite 
optical emissions (aurora) in the upper atmosphere. 

Both inside a cloud and in the plasma surrounding it, sec- 


ondary and higher-order curved charge layers are generated 
extending out into the ambient plasma. These layers have 
been discussed qualitatively in terms of the acceleration of 
ambient and cloud electrons in the vicinity of lower-order 
layers. However, the characteristics in terms of spatial sep- 
aration, expansion velocity into the ambient plasma, effects 
on cloud breakup, etc., have not been determined and will 
require longer-duration and more detailed simulations. 

A cloud that is pushed suddenly across a magnetic field 
(or a neutral gas cloud that is instantaneously ionized) will 
initially radiate close to the cloud electron and ion hybrid 
frequencies. If the ionization is slow compared to the par- 
ticle gyration period, particles will gyrate incoherently, and 
the radiation from the cloud will be of a smaller magnitude. 
Thus radiation at ion frequencies is more likely to be gen- 
erated through a gyrophase-coherent process than radiation 
at electron frequencies. Even when particles are random in 
phase, electromagnetic radiation will be generated at the gy- 
ration frequencies. Our results indicate that a cloud where 
the polarization field is effectively shorted will be halted 
quickly and the particles will gyrate as individual particles 
with almost their full initial perpendicular velocity. Recent 
observations in the CRIT 1 sounding rocket experiment of 
the CIV process [Brenniny et a/., 1991] included strong oscil- 
lations on the cloud ion cyclotron frequency and have been 
interpreted along the same lines as given above. 

In addition to the radiation from internal currents, cur- 
rents to the charge layers from the ambient plasma are ex- 
pected to radiate. In the Earth’s environment, the whistler 
mode should be generated through either Cherenkov reso- 
nance or cyclotron resonance. For the radiation to be signif- 
icant, the cloud dimension perpendicular to the velocity (y 
direction) must be large compared to a perpendicular wave- 
length. In this case radiation should be observed as in ar- 
tificial beam ejection experiments [Gurnett et a/., 1986] and 
naturally occurring beams in the auroral zone [James, 1976]. 
We note here that wave observations are potentially strong 
diagnostic tools for the probing of cloud dynamics provided 
the radiation aspects of clouds are well understood. The ad- 
vantage is that the observations can be made from remote 
locations and do not need to be in situ. 

Finally we point out the analogy of the electrodynam- 
ics of large conductors in space to cloud electrodynamics. 
For an insulated, perfect conductor with exposed conduct- 
ing end surfaces, the total electrostatic potential difference 
distributed between the two end surfaces and the ambient 
plasma is (? x 5) • L t where 0 is the velocity of the con- 
ductor and L is the vector from one end of the conductor 
to the other. This is the same expression as (4), which was 
used to calculate the electrostatic potential of the electron 
charge layer. The electromagnetic radiation and the current 
systems generated in the ambient plasma as a result of large 
orbiting structures and large plasma clouds have therefore 
much in common. The radiation from large space struc- 
tures have been studied by Donohue et ai [1991], current 
closure and radiation for large clouds have been treated an- 
alytically by Borovsky [1987], and current closure and wave 
generation for large space structures have been studied in 
laboratory plasma tanks by Stenzel and Urrutia [1990]. 
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Abstract We present here our first results of simulating 

the solar wind - magnetosphere interaction with a new three- 
dimensional electromagnetic particle code. Hitherto such globa 
simulations were done with NIHD codes while lower-dimensional 
particle or hybrid codes served to account for microscopic pro- 
cesses and such transport parameters as have to be introduced 
ad hoc in MHD. Our kinetic model attempts to combine the 
macroscopic and microscopic tasks. It relies only on the Maxwell 
curl equations and the Lorentz equation for particles, which 
are ideallv suited for computers. The preliminary results shown 
here are for an unmagnetized solar wind plasma streaming past 
a dipolar magnetic field. The results show the formation of 
a bow shock and a magnetotail, the penetration of energetic 
particles into cusp and radiation belt regions, and dawn-dusk 
asymmetries. 


I. Introduction 


T HE solar wind interaction with the earth’s magnetic 
fields sives rise to a number of important and intriguing 
phenomena, many of which are only partially understood. 
These include reconnection between the solar w’ind magnetic 
field at the day -side magnetopause and flux transfer events, 
the dras of the solar wind exerted on the magnetotail and 
associated instabilities at the magnetopause, plasma convec- 
tion in the magnetosphere/ionosphere, and the generation of 
field -aliened current systems. Reviews of these and other 
phenomena can be found in [l]-[4]. The range of physical 
processes involved in the solar wind -magnetosphere system is 
quite overwhelming, and consequently a wide array of methods 
have been used to study these, ranging from detailed studies 
of select phenomena with the assumption of, for instance, a 
specific field geometry or set of boundary conditions, to fully 
three-dimensional simulations with less spatial resolution but 
self-consistent particle or field geometries and with no (little) 
influence of the boundary conditions. 

The first three-dimensional (3-D), global magnetohydrody- 
namic (MHD) simulations of the solar w ind -magnetosphere 
svstem were reported in [5]-[7]. Since then, global 3-D 
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MHD simulations have been used with some success to 
study a range of processes including reconnection in the 
tail [8] -[9], reconnection at the day-side magnetopause [10], 
the dependence on the interplanetary magnetic field (IMF) 
orientation of magnetospheric convection and field-aligned 
currents [ 1 1 ]-[ 1 3], the solar wind - magnetosphere - ionosphere 
current -voltage relationship [14], the self-excitation of auroral 
arcs [15], and the reconnection voltage between the closed 
geomasnetic field and the IMF as a function of the IMF clock 
angle [16]. 

In MHD codes the microprocssses are represented by sta- 
tistical (macroscopic) constants such as diffusion coefficients, 
anomalous resistivity, viscosity, temperature, and the adiabatic 
constant. Recent studies have pointed out that results from 
MHD codes are sensitive to the choice of these parameters, as 
well as to the accuracy of the iteration scheme and have cast 
doubt on some results obtained in the past. Thus it has been 
found that in ideal MHD simulations using a high-precision 
code, no Kelvin -Helmholtz instability is generated along the 
magnetopause and there is no entry of mass, energy, and mo- 
mentum through the magnetopause, contrary to conventional 
knowledge [17]. 

To include particle dynamics explicitly, a number of simula- 
tions have been performed using hybrid codes (fluid electrons, 
particle ions) and full particle cedes in one and two dimen- 
sions. With such codes local problems have been investigated 
such as shocks formed in the solar wind as it encounters the 
magnetosheath [18]— [23] and the current layer at the day- 
side magnetopause [24]— [25]- Recently also global problems 
have been investigated in three dimensions with hybrid codes, 
namely the solar wind interaction with the day side of Venus 
[26] and the solar wand interaction with Mars and Venus [27]. 

With the model presented here, we intend to take the full 
step, namely the global simulation of the solar wind interaction 
with a planetary magnetic field using a particle code which 
contains the complete particle dynamics. As will be shown 
below, the advantage is that the model contains the complete 
physics. The price to be paid is that some scaling of plasma 
parameters must be done. 

II. A New 3-D EM Particle Simulation Model 

Our code is a successor to the TRISTAN code [28], [29]. 
Its new features are that (1) Poisson’s equation and Fourier 
transforms have been eliminated by updating the fields locally 
from the curl equations and depositing the particle currents 
according to charge -conserving formulas [30]; (2) radiating 
boundary conditions are applied to the fields using a first-order 
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Lindman approximation [31]; (3) filtering is done locally, (4) 
localization makes the code ideally suited to modem parallel 
machines which call for minimizing data paths; and (5) the 
code is in Fortran and is fully transportable: modest versions 
run on PC’s and on workstations. In the past, the TRISTAN 
code has successfully simulated large-scale space plasma phe- 
nomena such as the formation of systems of galaxies [3-]. 
The new version of the code has been applied to the study 
of the dvnamics of low - 3 plasma clouds [33] and a tutorial 
on the code will appear in the proceedings of the ISSS-4 
(International School of Space Simulation) held in Kyoto, 
Japan, in 1991 [34], 

For the simulation of solar wand- magnetosphere interaction 
the following boundary conditions were used for the particles: 
(1) Fresh particles representing the incoming solar wind (in 
our test run unmagnetized) are continuously injected across 
che yz plane at x = x min with a thermal velocity plus a buLk 
velocity in the -Fx direction; (2) thermal solar particle flux is 
also injected across the sides of our rectangular computation 
domain; and (3) escaping particles are arrested in a buffer 
zone and redistributed there more uniformly by making the 
zone conducting in order to simulate their escape to infinity. 
They are then written off. 

For the fields, boundary conditions were imposed just 
outside these zones: radiation is prevented from being reflected 
back inward, following Lindman’s ideas [31]. The lowest 
order Lindman approximation was found adequate: radiation 
at glancing angles w r as no problem. However, special attention 
was givenu> conditions on the edges of the computational box. 

In order to bring naturally disparate time and space scales 
closer together in this simulation of phenomena dominated by 
ion inertia and magnetic field interaction, the natural electron 
mass was raised to 1/16 of the ion mass and the velocity of 
light was lowered to twice the incoming solar wind velocity. 
This means that charge separation and anomalous resistivity 
phenomena are accounted for qualitatively but perhaps not 
with quantitative certainty. Likewise, radiation related phe- 
nomena (e.g. whistler modes) are covered qualitatively only. 


III. Results of First Tests 

While on a CRAY-2, with 100 M word core memory, a 
200 3 end simulation with some 10 million particles would be 
possible in principle, our tests were run on the CRAY-YMP 
at NCAR using a modest 105 by 55 by 55 grid and only 
200000 electron-ion pairs. Initially, these fill the entire box 
uniformly and drift with a velocity v D = 0.5c in the -i-x 
direction, representing the solar wind. The electron thermal 
velocity is v th = 0.02c while the magnetic field is initially 
zero. A circular current generating the dipole magnetic field 
is increased smoothly from 0 to a maximum value reached at 
time step 65 and kept constant at the value for the rest of the 
simulation. The center of the current loop is located at (40.5, 
27.5, 28) with the current in the xy plane and the axis in the z 
direction. The initial expansion of the magnetic field cavity is 
found to expel a large fraction ot the initial plasma. The solar 
wind density is about 0.7 electron-ion pairs per cell, the mass 
ratio is m l /m e = 16, and = 0.84. 


In Fig. 1 is shown the ion density in the center xz layer 
containing the dipole center at time step 785. (Data were 
recorded at 128 or 256 step intervals and while no fully steady 
state was approached, no major change was seen between 
steps 768 and 1024.) The ion density is color coded and the 
magnetic field component in the plane is shown with arrows 
at every third grid point. The magnitude of the field has been 
scaled in order to see the field direction for weak fields. Thus 
the length of the vectors is not a true representation of the 
field magnitude. The plasma is flowing through the simulation 
domain from left to right (low to high x values). In the process 
the dipolar field is compressed at the side facing the plasma 
wind and is extended to a long tail on the downwind side, just 
as the earth's magnetic field in the solar wind. Note, that in 
these simulations the current in the loop creating the dipole is 
such that the orientation of the dipole is opposite to that of the 
earth’s field. Or, the orientation of the plane of representation 
is different for what is normally used, with the North Pole 

down and the South Pole up. 

Some panicles have entered the cusp regions, and on the 
downwind side (night side), panicles have gathered on dosed 
field lines in what could be equivalent to the radiation belt. 
The panicles at the upwind magnetopause are lomning a bow 
shock. The temperature of these panicles is much elevated 
from the background temperature. Also seen upwind from the 
magnetopause is what could be a second shock or a foreshock. 
The magnitude of both shocks has a minimum at the subsolar 
point and is increasing to a maximum at some distance from 
the sun-eanh line. 

The current density associated with the electron and ion 
fluxes and the magnetic field are shown in Fig. 2. The current 
density component shown is perpendicular to the xz plane 
[J,). with current into the plane (the positive y direction) 
colored red and current out of the plane colored blue. Again, 
the values are for the center layer containing the center of the 
dipole. The current densities have been determined from the 
y component of the particle fluxes at each cell in the layer. 

The current system shown in Fig. 2 is equivalent to the 
Chapman-Ferraro current system of the earth. At the upwind, 
low - latitude magnetopause it runs in a direction w'hich cancels 
the original magnetic field on the upwind side, and doubles 
the field at the downwind side. At the cusp regions the current 
changes polarity with the magnetic field, as expecte-. The 
magnetopause in the tail also has negative J y . with a positive 
J y in the neutral sheath. The current in the neutral sheath is 
carried by the very few particles located here. 

The ion density and magnetic field in the yz plane are 
shown in Fig. 3. The values are those of the layer which 
contains the center of the dipole and the dipole axis. The 
direction of the solar w'ind is out of the plane toward the 
viewer. Note that the density is asymmetric, with more 10 ns 
located at the left side. This is caused by ions grad - B drifting 
toward this side as they move around the magnetic obstacle. 
The current density J z in the yz plane containing the dipole 
axis is shown in Fig. 4. The current is running in a belt around 
the dipole with a polarity that cancels fields on the outside, 
thereby creating a magnetopause. 
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Fig. 3. Ion density in the center y: layer containing the dipole center. 


Some key parameters are shown in Fig. 5 as functions 
of x. for center y and center z (sun -earth line). These are 
the ion density (top panel), the ion thermal pressure (second 
panel), the ion dynamic pressure (third panel), the magnetic 
pressure (fourth panel), and the ion Larmor radius (bottom 
panel). The magnetic field is small in the wind and in the 
tail, and rises to a peak at the location of the dipole center. 
Upwind, the increase is quite sudden and corresponds to the 
magnetopause. Upwind of the magnetopause, the density is 
hi eh and some oscillations in the density are present here. The 
nature of these oscillations has not yet been investigated. At 
the magnetopause the density drops, but recovers and reaches 
a local maximum at x = 39. This maximum is associated 
with cusp particles that have penetrated to the center of the 
dipole (there is no planet to stop them). At x = 46 a second 
maximum is reached which correspond to the radiation belt - 
tvpe particles mentioned earlier. Finally, the density increases 
slowly with distance in the tail. 

The thermal pressure increases dramatically in a region 
upw ind of the magnetopause. The location of this maximum 
corresponds to the location of the third of the three density 
maxima found upwind. Through the shock, the wind is slowed 
down, as seen by combining the plots of density and dynamic 
pressure. Thus the structures at the upwind side of the magne- 
topause bear much resemblance to the structures of the earth s 
bow shock. The ion Larmor radius has large variations w-ith 


location, as does the magnetic field. In the wind and in the 
tail (we are in the neutral sheath) the Larmor radius is of the 
order of 10-30, with peaks in places where the magnetic field 
becomes small. The magnetic fluctuations in the tail are caused 
by the tail waving like a flag. (In regions of rapidly varying 
low field strength, “Larmor radius" loses its meaning.) 


IV. Discussion 

The results presented in Figs. 1 through 5 show that even 
with the modest grid size of 10c by 55 by 55 cells, our a-D 
fully kinetic model is able to generate the complete magne- 
tosphere with its basic characteristics. This encourages us to 
continue our work using finer gnds (larger array sizes) and 
including such features as an IMF. a tilt of the magnetic axis, 
and an inner conducting surface simulating the ionosphere. 
Larger dimensions will allow finer resolution of fields. The 
number ot panicles should also be increased. Longer runs 
would allow’ us to enhance realism by using smaller m e /m t 
and larser c/vr>. On presently available single processor 
computers the grid resolution could be increased by a factor 
of 3 in all dimensions. Multiprocessor machines in the future 
will allow even greater improvement, since the code allows 
simultaneous updating of many cells. 

Some important issues that with advantage can be addressed 
with a kinetic code are as follows: 
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1) In MHD simulations reconnection phenomena are some- 
times studied using artificial numerical dissipation to 
drive the process [16]. Furthermore, artificially high 
diffusion coefficients are often used in order to achieve 
a stable solution [11]- With a full kinetic code such 
problems may be addressed directly. 

2) The motions of individual panicles are tracked. In many 
regions of the magnetosphere, such as the magnetopause 
and the tail, gradients in the magnetic field are of 
such magnitude that the magnetic field changes sig- 
nificantly during a panicle gyration. This gives rise 
to nonisotropic, non-Maxwellian distribution functions 
which drive the physics of those regions. The recent 
research into chaotic panicle orbits in the magnetotail 
[35*. [361 is an example. The complex panicle motion 
ci\es rise to such phenomena as collisionless electron 
viscosity [37]— [39] and collisionless conductivity [40], 
Furthermore, a panicle code calculated the panicle dy- 
namics in self-consistent rather than model magnetic 
fields. 

3) The complete wave dynamics from low-frequency mag- 
netohydrodvnamic modes (Kelvin — Helmholtz instabili- 
ties at the magnetopause and ULF pulsations) to electron 
waves (Langmuir oscillations) are included. An example 
where electron kinetics is important is illustrated by 
recent studies of the dynamics of the magnetopause 
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inter. Positive currents f.ow out of the plane and negative into the plane. 

current layer [25] and the dynamics of expanding plasma 
clouds [41]. It is thought that a generalized lower hybrid 
drift instability is generated, which in the case of the 
magnetopause has the effect of broadening the current 
layer and causing filamentation of the current. 

4) Any complicated magnetic field configuration can be 
simulated by a simple modification of the source current. 
This makes it straightforward to simulate not only other 
planetary' magnetospheres such as that of Uranus, but 
also the effect of the diurnal rotation of the earth's dipole 
axis with respect to the direction to the sun. 

Questions that need to be addressed in the future are those 
of scaling from the microworld of the simulation to the real 
macroworld, in particular, the scaling of the ion to electron 
mass ratio, the Larmcr radii, and the Debye length to the 
earth’s radius or other characteristic dimensions. An example 
is that in our simulation the Debye length is a significant 
fraction of the earth’s radius or the thickness of the bow 
shock. The earth's radius itself can be determined by the 
distance from the dipole center to the upwind magnetopause. 
This distance is typically equal to ten earth radii. In the test 
simulation described earlier an earth radius is of the ordei of 
one grid spacing. 

However, some form of scaling is usually needed in particle 
simulations, even in one and two dimensions, and still such 
simulations are able to reveal much of the physics behind 
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Fig. 5. Key parameters as a function of .r at center y and center r (along 
the sun -earth line). From top to bottom arc: density, the thermal pressure, the 
dynamic pressure, the magnetic pressure, and the ion Larmor radius. 


natural phenomena. Similar problems are also encountered, 
and to some degree overcome, in laboratory plasma simu- 
lations of macroscale space plasma phenomena. Examples 
include the study of magnetic held line reconnection [42], 
[43] and the interaction of a magnetized plasma flow with 
a dipole magnetic field [44]. We propose, therefore, that a 
kinetic particle model will be a useful tool for the study of 
the earth's and other planetary magnetospheres in the solar 
wind, by supplementing present models that do not contain 
the complete panicle dynamics. 
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Abstract We have restudied electron beam driven whistler 
waves with a 3-D electromagnetic particle code. In the ini- 
tialisation of the beam-plasma system, U quiet start” conditions 
were approached by including the poloidal magnetic field due 
to the current carried by beam electrons streaming along a 
background magnetic field. The simulation show electromag- 
netic whistler wave emissions and electrostatic beam modes like 
those observed in the Spacelab 2 electron beam experiment. It 
has been suggested in the past that the spatial bunching of 
beam electrons associated with the beam mode may directly 
generate whistler waves. However, the simulation results indi- 
cate several inconsistencies with this picture: (1) the parallel 
(to the background magnetic field) wavelength of the whistler 
wave is longer than that of the beam instability, (2) the paral- 
lel phase velocity of the whistler wave is smaller than that of 
the beam mode, and (3) whistler waves continue to be gener- 
ated even after the beam mode space charge modulation looses 
its coherence. The complex structure of the whistler waves in 
the vicinity of the beam suggest that the transverse velocity 
(gyration) of the beam electrons in the combined poloidal and 
background magnetic fields is involved in the generation of the 
whistler waves. 


1. Introduction 

During the Spacelab 2 flight, launched on July 29, 1985, a 
spacecraft called the Plasma Diagnostics Package (PDP) per- 
formed a fly-around of the shuttle at distance of up to 300 me- 
ters while an electron beam was being ejected from the shuttle 
in a steady (d.c.) mode [Gumett et a/., 1986]. The plasma 
waves observed as the PDP approached and crossed through 
the beam took the appearance of a well-defined, funnel-shaped 
feature when displayed on frequency-time spectrograms. The 
waves were observed over a broad frequency range extending 
up to the electron cyclotron frequency (/ c ft? 1 MHz) with the 
center of the funnel coinciding with the passage through the 
beam. The mode was identified as the whistler mode, propa- 
gating with the wave normal close to the resonance cone. At 
even higher frequencies (3.1 MHz), intense narrowband emis- 
sions were observed inside the beam. These emissions are be- 
lieved to be either at the local electron plasma frequency, f p 

or the upper hybrid resonance frequency, Juhh — ( fc + fp )^ 2 
[Gumett et ai , 1986]. 

The generation of whistler waves by natural and artificially 
injected electron beams has been studied extencively in the past 
by means of space experiments, analytical theories, and particle 
simulations. A review of recent work is given by Neubert and 
Banks [1992]. However, the relationship between the electro- 
static beam mode (bunched beam electrons) and the whistler 
waves is not yet clearly understood. The early work by Bell 
[1968] proposed that Langmuir turbulence in natural auroral 
electron beams provided a means for quasi-coherent Cherenkov 
radiation needed to explain the relatively high level of whistler 
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wave radiation observed in the auroral regions. In the Spacelab- 
2 experiment it was also found that some mechanism creating 
quasi coherence was needed to explain the observed signal lev- 
els. Employing 1-dimensional electrostatic particle simulations, 
Farrell et ai [1989] found that Langmuir turbulence in the 
beam created a spatial modulation of the current structure, 
and therefore enhanced the the ability of the beam to radiate. 
The estimates of the radiation levels based on the simulations 
were found to be in good agreement with observed levels of 
radiation. Electromagnetic simulations in 2 dimensions con- 
firmed this picture [Pritchett et ai, 1989]. Here it was found 
that even with no initial net current (background electrons drift 
opposite to beam electrons thereby providing a return current 
to a spacecraft), instabilities are generated inside the beam 
creating a current structure which allow for the radiation of 
electromagnetic waves. 

The outstanding question is really the generation mechanism 
of the radiating current structure, its dynamics, and how it de- 
pends on ambient and beam plasma conditions. Several studies 
have appeared which imply that the simple picture of an elec- 
tron beam, space charge modulated by Langmuir turbulence, 
is not the means through which whistler radiation is generated, 
although it may be an important “seed” mechanism to get the 
process started. These works include Goerke et ai [1990] who, 
in a sounding rocket experiment, found no evidence to sup- 
port the hypothesis that electrostatic bunching at the plasma 
frequency was responsible for the enhancement of the whistler 
mode radiation. Furtheremore, simulation results have pointed 
out that the wavelength of the electrostatic beam turbulence 
differs from the wave length of the whistler wave [Omura and 
Matsumoto, 1988a]. 

In this work we attempt to establish more insight into the 
generation mechanism of whistler waves using a newly devel- 
oped 3 -dimensional, electromagnetic, and relativistic particle 
code [Buneman, 1992; Buneman et al , 1992; Neubert et ai, 
1992]. It extends the 3-dimensional simulations of Nishikawa 
et ai [1989] by using a larger system, by including all compo- 
nents of the magnetic fields as well as relativistic effects, and by 
including a self-consistent poloidal magnetic field as an initial 
condition. 


2. Simulation Results 

The system size used for the simulations is L x = L y = 
85A,L* = 160A, where L x , L y , and L z are the lengths of 
the system in the three dimensions and A(= 1) is the grid size. 
Periodic boundary conditions are used along the ambient mag- 
netic field (z-direction), while radiating boundary conditions 
[Lindman, 1975] are used in the x- and y-directions. 

The electron beam is initially in a column with radius, r e 5 = 
4.47A, and the axis at the center of the domain (x = y = 43) 
oriented in the z-direction. The plasma density within the 
beam column is the same as outside and the beam is initially 



3 


charge neutral. The location of beam electrons within the col- 
umn are randomly chosen and is about half of the total elec- 
tron density. The total number of beam electrons is 4899. 
The beam drift velocity, = 3.5i; e t, where v ct is the elec- 
tron thermal velocity. The beam electron temperature, = 
0.097; , where T e is the temperature of the background elec- 
trons. In this study the ion drift velocity, V x is zero. To approx- 
imately simulate ionospheric conditions we chose the following 
parameters: m,/m c = 64, T e JT{ = 1, Qe/vpe = 0.4, c/v et = 
10.67, ujpcAt = 0.10, 0 = 0.111, X Dc = 0.469A, p c = 1.17, 
Pi = 9.38, c/lj pc = 5.00. 

As a new feature, quiet start conditions are achieved by 
including an initial poloidai magnetic field (IPMF) consistent 
with the current J z carried by the beam electrons at t = 0. 
The IPMF is calculated from the curl of the vector potential, 
A, where A = (0,0, A z ) and A z is given by, 


A z (i,j ) = A zo \og(r 2 /r 2 eb )\r > r ei> , (1) 


A z (i,j) = A„(log(1.0) + r 2 /r 2 eb - 1.0); r < r eb , (2) 

where A zo = —J z /(4irc) ) r is the radius from the beam column 
axis, and J z is the current in the column. 

As mentioned earlier, Pritchett et al [1989] also simulated 
a quiet start ( J z = 0). A quiet start avoids a transient stage 
where a strong return current is induced as reported by Omura 
and Matsumoto [1988a, b]. 

In Figure la, B x , B y (with the IPMF subtracted) is shown in 
the center plane perpendicular to the beam axis at = 25.6. 
The length of the vectors have been scaled to enhance weaker 
signals in order to better illustrate the field structure. A com- 
plicated wave pattern resembling a honeycomb structure is seen 
rather than simple concentric wave fronts. This structure oscil- 
lates in time (not shown) constituting an outward propagating 
electromagnetic wave. 

Figure 1 b shows the magnetic field components B X) B Z in 
the center (x,z) plane. The ambient magnetic field B z0 has 
been subtracted. Waves are clearly created near the beam and 
propagating out in a complicated manner. If the parallel cur- 
rent is uniform, as for the initial conditions, the values of B x , B z 
should be zero in this plane. The fact that these components 
take on finite values in the center of the beam leads to the idea 
that a spiral structure of the current is produced in addition 
to the current along the z-direction, as will be discussed below. 
Obviously, the B z - component is created by the electron mo- 
tions perpendicular to the ambient magnetic field, which also 
creates transverse electric fields E X) E y . This is essential to 
whistler waves. 

The beam electrons in phase space z — v z is shown in Figure 
2a. The electrons are bunched, accelerated, and decelerated, 
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fundamentally due to the excited electrostatic waves. The back- 
ground electrons, initially in the column, also show bunching 
and acceleration in and around the beam column (not shown). 

As a result, a perturbed current structure is generated which 
acts like an antenna and emits whistler waves as suggested by 
[Farrell et a/., 1989; Pritchett et a/., 1989; Wong and Lin , 1990]. 

The ./^-component of the current is shown in Figure 2b. How- 
ever, new to our study is that altough the beam initially was 
directed along the ambient magnetic field with no transverse 
component, transverse motion develop which are es- 

sential to the whistler wave generation. 

The time evolution of the charge density and the magnetic 
field component B x along the beam axis is illustrated in Figure 
3. The charge density shows sporadic waves structures prop- 
agating in the beam direction (towards larger z). The phase 
velocity is of the order of 2.6 v e t for 12 < u p€ t <21. The av- 
eraged drift velocity of the beam electrons varies at the same 
time from 3.26 v e t to 2.56t> e t. The phase velocity of the wave is 
less than but of the order of the averaged drift velocity. There- 
fore, this wave is excited by the beam electrons and is identified 
as a beam mode. On the other hand, the magnetic field com- 
ponent B x , which represent a whistler wave, shows a longer 
wavelength and a more clear coherent structure. The phase ve- 
locity is about 2.22v ct for the same time period, which is slower 
than for the beam mode. 

The space charge modulation at the earlier time (u > pe t < 

10) does not excite whistler waves. From this fact we may 
draw the idea that the bunching of beam electrons becomes 
the seed of the excitation of whistler waves, and only when 
the electrons are affected by the poloidal magnetic field, which 
create a transverse velocity, are whistler waves excited. This 
process was first suggested by Bell and Buneman [1964]. 

Although difficult to conclude from Figure 3a, a closer in- 
spection of the simulation results indicate that the space charge 
oscillations later get randomized while whistler waves continue 
to be generated. This means that the whistler wave and the 
beam mode are well separated, a conclusion reached also by 
Omura and Matsumoto [1988a] and Goerke et a/., [1990], as 
mentioned earlier. It is also consistent with the idea that 
the self-quenching of the electrostatic instability by creation 
of an enhanced (turbulent) velocity spread may still allow the 
whistler mode to take over eventually [Bell and Buneman , 1964].| 

It suggests that some current structure still exist in the 3-dimensiona% 
system (helically). 

In Figure 4 the accumulated spectra along the x direction is 
plotted for the magnetic field component B x . The color scale 
is logarithmic in wave intensity. High-intensity waves are seen 
below the hybrid frequency within the electron beam, and lower 
intensity whistler waves below the electron cyclotron frequency 
propagateing outside the beam. 

The dispersion of the electromagnetic waves just outside of 
the beam column is shown in Figure 5. The plot has been 
created by performing a 2-dimensional Fourier transformation 
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(one space and one time) of B x along the ^-direction for the first 
512 time steps. The whistler branch extents from the origin to- 
ward the horizontal line at the electron cyclotron frequency Q c 
with additional, lower amplitude radiation generated at higher 
frequencies along a line with a slope corresponding to the ve- 
locity of light. 


3. Discussion 

We have attempted to establish a clearer picture of the 
generation mechanism of whistler waves. The simulation re- 
sults show an electrostatic mode contained within the beam 
collumn and whistler emissions propagating from the beam 
into the ambient plasma. These results agree qualitatively 
with observations obtained by PDP during the Spacelab 2 mis- 
sion. As expected, the simulation results have many similarities 
with those of previous 1- and 2-dimensional simulations. The 
initial homogenous parallel current is perturbed through the 
beam mode, developing spatial inhomogeneities, which stim- 
ulate whistler radiation [Farrell et ai , 1989; Pritchett et a/., 
1989]. On the other hand, we found that the beam mode is 
well separated from the whistler mode as also pointed out by 
2-dimensional simulations [Omura and Matsumoto , 1988a] and 
sounding rocket experiments Goerke et ai [1990]. 

Because the beam mode must be instrumental in the pertur- 
bation of the beam current, we are faced with the challenge of 
unifying seemingly contradictory evidence. We suggest a solu- 
tion along the lines first discussed by Bell and Buneman [1964]. 
Here it was found that a spread in the beam electron velocity 
in the transverse direction leads to an instability of the stream- 
plasma system in the whistler mode, with the transverse elec- 
tron gyrations serving as an energy source. As shown in Figure 
lb, the magnetic componenet B z is produced by the transverse 
electron motion which is consistent with this finding. 

While the electrostatic modes are saturated and quenched 
due to the slowdown and heating of beam electrons (around 
u pc t = 25), whistler waves are still being excited. At this 
stage the bunching of beam electrons has disappeared in the 
z — v z phase space. Phase mixing through longitudinal velocity 
spread can suppress the electrostatic instability without sup- 
pressing whistler growth as pointed out by Bell and Buneman 
[1964]. Therefore, a 3-dimensional spiral bunching structure of 
the beam (background) electrons may be responsible for this 
prolonged excitation of the whistler mode. Further discussions 
will be reported in a separate publication. 
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Figure Caption 

Fig. 1. The fields in the 2-dimensional planes with the initial 
poloidal magnetic field. 

(a) B x - By in the x — y plane ( z = 80) at u pe t = 25 

(b) B x - B z in the x — z plane (y = 43) at u p€ t = 25.6 

Fig. 2. Phase space distributions (z — v x of beam and the total 
current by the electrons initially in the column. 

(a) The beam electrons at u) pe t = 25 

(b) The total current by the beam and background electrons 
initially in the column at u pe t = 25 

Fig. 3. The time evolutions of the wave patterns along the z 
direction at x = 43, y = 43. 

(a) the charge density 

(b) the the magnetic field component ( B x ) 

Fig. 4. The spectrogram of the magnetic field componenet 
( B x ) along the x direction at y = 43, z = 80 for the period 
0 < < 51.2. 

Fig. 5. Wave intensity as function of u and k z along 22 < z < 
150 at x = 43, y = 20. 
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